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ABSTRACT 

We present the analysis of a kinematic data set of stars in the globular cluster NGC 2419, taken 
with the DEIMOS spectrograph at the Keck II telescope. Combined with a reanalysis of deep Hubble 
Space Telescope and Subaru Telescope imaging data, which provide an accurate luminosity profile of 
the cluster, we investigate the validity of a large set of dynamical models of the system, which are 
checked for stability via N-body simulations. We find that isotropic models in either Newtonian or 
Modified Newtonian Dynamics (MOND) are ruled out with extremely high confidence. However, a 
simple Michie model in Newtonian gravity with anisotropic velocity dispersion provides an excellent 
representation of the luminosity profile and kinematics of the cluster. The anisotropy profiles of these 
models ensure an isotropic center to the cluster, which progresses to extreme radial anisotropy towards 
the outskirts. In contrast, with MOND we find that Michie models that reproduce the luminosity 
profile either over-predict the velocity dispersion on the outskirts of the cluster if the mass to light ratio 
is kept at astrophysically-motivated values, or else they under-predict the central velocity dispersion 
if the mass to light ratio is taken to be very small. We find that the best Michie model in MOND 
is a factor of ^ 10'' less likely than the Newtonian model that best fits the system. A likelihood 
ratio of 350 is found when we investigate more general models by solving the Jeans equation with 
a Markov-Chain Monte Carlo scheme. We verified with N-body simulations that these results are 
not significantly different when the MOND external field effect is accounted for. If the assumptions 
that the cluster is in dynamical equilibrium, spherical, not on a peculiar orbit, and possesses a single 
dynamical tracer population of constant M/L are correct, we conclude that the present observations 
provide a very severe challenge for MOND. 

Subject headings: gravitation — globular clusters: individual (NGC 2419) — stellar dynamics 



1. INTRODUCTION 

Are accelerations in the weak regime well described 
by the Newtonian approximation to General Relativity, 
or is a modification to this the ory, such as th e Modified 
Newtonian Dynamics (MOND; jMilgrom 11983), required? 
MOND was proposed 30 years ago as an alternative to 
dark matter to explain the approximately flat rotation 
curves in the outer parts of disk galaxies. In this the- 
ory the gravitational acceleration deviates significantly 
from that predicted by Newtonian gravity only in the 
vicinity or below a value ao(^ 1.2 x lO~^cms~^). One 
implementation of this idea, due to [Bekenstein fc MlT] 
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grom ( 1984 1 , is to modify Poisson's equation in Newto- 
nian gravity 'S/'^ipN = AirGp to: 



(^)v. 



47rGp, 



(1) 



where <f is the MOND potential, and /i is an interpolating 
function that changes between fi{x) = 1 in the Newto- 
nian regime, and ^{x) = a: in the deep MOND limit of 
low accelerations. For a finite-mass isolated system the 
boundary conditions of equation |l]) are \^ip\ — > for 
\r\ — ;> cxD, where f is the position vector. 

While MOND was initially developed to explain the 
rotation c urves of ga laxies and the Tu Uy-Fisher relation 
(see, e.g., [TullyfeFisher 19771 pcGaugh 2005) , it was 
not evident that it would hold up to closer scrutiny in 
other astrophysical systems. Indeed, many observational 
challenges to MOND have been presented over the years. 
The theory has remained surprisingly resilient to this ob- 
servational assault (Sanders & McGaugh 2002), though 
in some cases there are dithculties in reconciling the ob- 
servations with the predictions of MOND. These include 
problems with the growth o f cosmological perturbations 
(Dodelson & Liguori 2006), the offset bet ween lensing 



mass and baryons in th e Bullet Cluste r ( Clowe et al . 
'2006'), Solar System tests ('Milgrom'1983'; Sereno & Jetzer 
2006 ) , dynamical fricti on in dwarf galaxies ( ^lotti & Bin- 
ney 2004;s Sanchez- Salcedo et al.l2006[|Nipo^e t al. 2003f 
I Angus fc Diaferio||2009 ), and the kinematics and den- 



sity profile of satellite galaxies (Klypin & Prada 2009), 
to list a few. However the main-stream explanations of 
these systems that involve dark matter are not them- 
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selves without problems, especially at the galactic scale 
where the discrepancy between the predictions of stan- 
dard Cold Dark Matter theory (ACDM) and observations 
appears to be greatest. The issue of the missing satel- 
lite galaxies, the possible non-existence of dark matter 
cusps and the high angular momentum of galactic disks 
are among the more famous problems that pose difficul- 



ties to ACDM (Binney 2004 



Primack||2009[ ), ahhough it 



is possible that the discrepancies are due to our incom- 
plete understanding of the physics of baryonic matter. 
It is therefore very important to challenge these theories 
of gravity in the most stringent ways we have access to, 
using well-understood astrophysical systems, that can be 
measured and probed with great accuracy. 



It has been appreciated for a number of years ( Scarpa 
et al.|2003[ |Haghi et al.|2009[ ) that the dynamics ot glob 



ular clusters may give a means to assess the fundamental 
question above, since the accelerations that stars expe- 
rience on the periphery of these systems are exceedingly 
low. Indeed, if one were to consider a globular cluster in 
isolation, one would expect to see a drop in the acceler- 
ations that stars experience down to values approaching 
zero on the periphery of the cluster. In reality, however, 
these systems reside inside a host galaxy, so in general 
they cannot be treated as isolated. In Newtonian dynam- 
ics, provided that the gravitational field of the galaxy can 
be considered uniform across the extent of the cluster 
(i.e. if tidal effects are negligible) , the internal dynamics 
of the cluster are the same as if the cluster were isolated. 
This is not the case in MOND, in which, due to the non- 
linearity of the theory, the internal dynamics of a stellar 
system is influenced even by a uniform external field (in 
other words, th e strong equivalence principle is not valid; 
Milgrom 1983 ). In the considered case of MOND as mod- 
ihed gravity, this external field effect (EFE) can be ac- 
counted for in the boundary conditions of Equation ([T]) , 

which are V(/3 — >■ V^pext for |r| — > oo in the presence 

of a uniform ex ternal field g^xt — —^^ext (Bekenstein & 
Milgromjl984 ). An important consequence ot the El^'E is 



that the internal dynamics of a MOND system with even 
low internal accelerations becomes almost Newtonian in 
the presence of a sufficiently strong external field. It is 
then clear that for a globular cluster to be useful in the 
attempt to distinguish between Newtonian gravity and 
MOND, not only must its internal accelerations be low, 
but also the external galactic field it experiences must be 
low, compared to cq. 

Figure [iji illustrates the external gravitational acceler- 
ation due to the Milky Way that members of the Galactic 
globular cluster population experience. For the majority 
of the clusters, the external field accelerations are sub- 
stantially above the MOND characteristic acceleration 
ttQ. Note also that the Galactocentric distances shown in 
Figure [T^ refer to the present-day orbital phase, so that 
even some of the more distant clusters will pass closer 
to the Milky Way center during their orbits. It is only 
the most distant clusters, well in the realm of the distant 
halo, that can be considered to be sufficiently isolated to 
probe the low acceleration regime. Also relevant is the 
fraction of the cluster stars that lie within the MOND 
regime: this fraction is substantial only for the most ex- 
tended clusters (Figure [Tb). 

NGC 2419 is a particularly promising globular cluster 
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Figure 1. Top panel: The expected external accelerations that 
Galactic globular clusters experience as a function of their radial 
distance from the ce nter of the Milky Way. Th e realistic Galac- 
tic mass model from [Dehnen &: Binney| l|1998[ l is assumed, and 
the actual three-dimensional positions of the clusters are used. 
Though the acceleration is computed for a Newtonian model of the 
Milky Way with dark matter, we expect similar values in MOND 
if MOND has to reproduce the rotation curve of the Milky Way. 
Those globular clusters that experience an acceleration below ao 
are circled red, and the subset of these clusters that are expected to 
contain fewer than 50 RGB stars are further marked with crosses. 
The crossed-out clusters are poor targets since they contain very 
fe w stars amenable to spectroscopic follow-up (see the discussion 
in |Haghi et al.|2009[| . Clearly, NGC 2419 is the only rich cluster in 
the outer Galaxy. (We have shifted the two clusters on either side 
of NGC 2419 — Eridanus and Palomar 3 — by 5 kpc to avoid a 
confusing overlap at Rqc ^ 95 kpc). Bottom panel: the expected 
internal acceleration at the half-light radius, assuming Newtonian 
gravity, for the same sample of globular clusters, as derived from 
the compilation by Harris (1996). Clearly NGC 2419, the target of 
this study, is minimally affected by the external Galactic accelera- 
tion field, and is also sufficiently extended so as to have a significant 
amount of its member stars in the low acceleration regime. The 
MOND ao level is indicated with a dotted line in both panels. 

to study in this context. First, as indicated in the top 
panel of Figure [TJ it is one of the most distant objects 
among the known population of Milky Way clusters, so 
it experiences a low external force. Indeed, the accelera- 
tion due to the Gala xy at that location , accor ding to the 
Milky Way model of Dehnen & Binney ( 1998 1, is well be- 
low the critical MOND acceleration value [gext ~ 0.1 ao). 
Furthermore, NGC 2419 is a very extended cluster, which 
implies that the accelerations on its constituents are, on 
average, very low. Thus is illustrated in the lower panel 
of Figure [l] where the internal acceleration is given as 
a function of the half light radius. Clearly the expected 
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Table 1 

Properties of the globular cluster NGC 2419. 
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Figure 2. The Newtonian internal accelerations in the best-fitting 
Michie model of NGC 2419 found in this study (see §4), as a func- 
tion of distance from the cluster center. The core radius rc, the 
half-light radius r^, and the tidal radius rt are marked. With our 
improved parameter constraints, it is evident that even at the clus- 
ter half-light radius, the acceleration that stars experience is below 
the ag value. 

accelerations drop rapidly further out in the cluster, and 
as we show in Figure [2[ slightly beyond the half mass 
radius in NGC 2419 fall well below gq. The high lumi- 
nosity of this cluster (which is the fourth most luminous 
in the Milky Way, after w-Cen, M54 and NGC 6388) 
also means that there are large numbers of target stars 
for study, even in the shortest evolutionary phases, such 
as the the upper Red Giant Branch (RGB), and also 
ensures that the velocity dispersion is large enough to 
resolve easily with current instrumentation. These prop- 
erties make NGC 2419 a muc h more favourable subject 
than nearby massive clusters ( Lane et al.|2010l), or dis- 
tant low-mass clusters like Palomar 14 ( jjordi et al.|2(J(J9 
Gentile et al.|2010 ) that have previously been scrutmizec 
to test MUND. 

Throughout this study, unless otherwise indicated, we 
adop t the distance measuremeii t of (m — M)o — 19.71 ± 
0.08 (Di Criscienzo et al. 2011); this implies a distance 
of 87!5kpc and an angular scale of 25.452 pc per ar- 
cmin. We assume the center of the cluster is located 
at 07''38™08^47 -t- 38°52'55".0 (iDalessandro et al.||2008l 



^08^47 -t- 38°52'55".0 ( |Dalessandro et al „ ^ 

hereafter D08), which is towards the Galactic anti-center, 
where the reddening is measured to be E(Ti—y) = 
0.08±0.01 ( [Ripepi et al.|2007|p^riscienzo et al.|2011 1 
from RR Lyrae lightcurves. 

2. OBSERVATIONS 

The analysis in this article is based upon several photo- 
metric and spectroscopic data sets. The main photomet- 
ric data set, which serves to determine the stellar surface 
brightness profile of the cluster, is derived from images 
taken with the Advanced Camera for Surveys (ACS) on 
board the Hubble Space Telescope and much wider-field 
images taken with the Suprime-Cam camera on the Sub- 
aru Telescope, as derived by D08. During the course 
of the analysis presented below, we realized that D08 
had overestimated the uncertainties on the density pro- 
file, making the model comparisons with the data highly 
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0.08 ± 0.01 mag 
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87.5 ± 3.3 kpc 
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1.90 ± 0.16 Mq/Lq 
9.12 X 10^ 
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surprising in the sense that they were "too good to be 
true" . To address this problem we used their photomet- 
ric catalog to re-determine the light profile of the clus- 
ter. We define a sub-sample of 7849 stars, by selecting 
from a color-magnitude box that incorporates the up- 
per main sequence, the RGB and blue horizontal branch 
(BHB) populations, and limited the sample to V = 23.5 
(as argued by D08) to avoid incompleteness, and radial 
variations thereof. 

Star counts were obtained in 15 radial annuli (or por- 
tions of radial annuli) within r < 612", the largest circle 
that is fully included in the Suprime-Cam field. As can 
be appreciated from Figure 2 of D08, although the Sub- 
aru data cover the cluster well, they do not extend suffi- 
ciently far away from the cluster so as to provide a back- 
ground comparison region that is unambiguously free of 
cluster member stars. Accurate background subtraction 
is a crucial issue for our study, so to explore better the 
transition between the cluster and the background popu- 
lation^, we obtained archival images of the ROSAT clus- 
ter RXCJ0736. 4-1-3925 taken with the MegaCam Camera 
on the Canada- France- Hawaii telescope (CFHT). With 
exposures of 2160 sec in g and 1920 sec in the r-band, 
these images cover a 1° x 1° field, with NGC 2419 lo- 
cated on the southern edge of the mosaic. The MegaCam 
images were processed and reduced to photometric cata- 
logs in an identical way to that described in Ibat a et al.| 
(20071, u sing the Cambridge Astronomical Survey Unit 
pipeline (Irwin & Lewis 2001). By cross- identifying the 
sources in the Cl^'Hl' and Subaru catalogs, we devised a 
color-magnitude selection in the g and r bands that was 
as close as possible to that made in the V and I bands 
from the D08 sample. Figure [3] shows the star-counts 
profile of the point sources in the CFHT image that lie 
within the defined stellar locus. Clearly, beyond 10' and 
up to 1° from the cluster center the star counts are con- 
sistent with a flat background. This justifies our choice 
of using the area at r > 612" in the Suprime-Cam/ACS 
survey (our main photometric catalog) to estimate the 
background count level (0.00013 ± 0.00002 arcsec'^). 

The absolute photometric zero-point of the Suprime- 
Cam/ACS profile was defined by matching the data at 

^ The "background" population in reality is of course mostly due 
to foreground Galactic stars. 
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Figure 3. The star-count profile in tiie CFHT field used to probe 
the background profile at large distance from the cluster. The 
stellar selection was chosen to be as close as possible to that 
of the main Suprime-Cam/ACS survey, and indeed the weighted 
mean background level of the CFHT counts between 20' and 60' is 
0.00014 arcsec"^ (dashed line), very similar to the Suprime-Cam 
average of 0.00013 ± 0.00002 arcsec-2 (dotted fine) at r > 612". 
The CFHT data shows that there is no significant contribution 
from the cluster at radii beyond ~ 10'. 

r < 16"5 with t he 6 innermost points of the profile by 
Bellazzini (2007 hereafter B07) that are based on aper- 
ture photometry of the ACS image. The gaps between 
detectors and the most obvious holes due to saturated 
stars were taken into account in computing the area of 
the annuli. The resulting profile is tabulated in Table [2j 
apart from the improved uncertainty estimates, this pro- 
file is in agreement with that of DOS. Note, however, that 
the non-azimuthally symmetric geometry of the missing 
areas (see Figures 1 & 2 of DOS) gives rise to slight ra- 
dial biasses in the star-counts profile, which must be ac- 
counted for in any fitting or modeling of the cluster by 
the use of appropriate window functions^. 

Targets for follow-up spectroscopy were selected from 
an earlier compilation of photometric data constructed 
by B07, which is a combination of photometry measured 
from the same ACS images, as well as from ground- 
based images which cover a wider field of view. The 
ground-base d phot ometry includes the WIYN study by 
Saha et al.| ([gOOSi Sloan Dig i tal Sk y Survey (SDSS) 
data [Adelman-McCarthy et al. (2006) as well as images 
retrieved from the archive ot the i'eiescopio Nazionale 
Galileo (TNG). 

The main spectroscopic observations presented here 
were t aken with the DEI MOS multi-object slit spectro- 
graph (Faber et al. 20031 mounted on the Keck II tele- 
scope at Mauna Kea, Hawaii on 2008 October 27th. We 
employed the high- resolution 1200 line/mm grating cou- 
pled with the OG550 order-blocking filter to give access 

o 

to the spectral range ~ 6500 - 9500 A at a dispersion 

of 0.33 A/pixel (i.e. 12kms-Vpixel at S500 A). The 
spectrograph has a scale of 0"1185/pixel and an anamor- 
phic magnification of 0.7. The spectrograph camera has 

* The window functions and color-magnitude selection functions 
are available from the authors, on request. 



two rows of detectors 'blue' and 'red'; in order to main- 
tain the Ca II triplet lines on the red side of the camera 
to avoid the inter-CCD gap, we set the central wave- 

o 

length to 8000 A. Narrow slits of 0"7 were milled to op- 
timize the spectral resolution, which was determined to 

o 

be ~ 1.2 A, based upon measurements of the width of arc 
lamp lines. Target stars were selected from the B07 cat- 
alogue with I < 20 and within a narrow color-magnitude 
region around the cluster Red Giant Branch (RGB), as 
shown in Figure |4]3. Two masks milled with slitlets for 
85 and 77 targets, were observed with exposure times of 
3 X 1200 s and 3 x 1000 s, respectively (we will henceforth 
refer to these as masks 1 and 2). The weather was clear 
during these exposures and the seeing was approximately 
0'!7, fortuitously matching the slit widths. 

In order to convey the quality of the data, we display in 
Figure [5] the spectra (and the corresponding noise spec- 
tra) for three stars that are representative of the bright- 
est, of intermediate and of the faintest stars in the clus- 
ter sample. The signal to noise of the spectra are clearly 
high, exceeding S/N = 10 per pixel for the faintest stars 
in the sample and reaching as high as S/N ^ 70 for the 
very brightest stars. 

A new spectral reduction package was constructed, 
partly for this project, in an attempt to improve object 
extractions, wavelength calibration and sky subtraction, 
as well as to provide more reliable uncertainty estimates. 
The procedure will be presented in detail in a forthcom- 
ing contribution (Ibata et al., in preparation); a brief 
summary is sufficient for the purpose of the present work. 
The aim of the software is to avoid any resampling of the 
observed two-dimensional spectroscopic frames: each ob- 
served pixel and its corresponding uncertainty is kept 
all the way to the final calibrated product, at which 
point the spectra may be trivially collapsed into one- 
dimensional form for visualization purposes, if necessary. 
The reduction steps are as follows: 

• identify cosmic rays in the images if multiple expo- 
sures are available. 

• correct for scattered light. 

• perform a flat-field correcting for pixel-to-pixel 
variations. 

• perform an illumination and the slit function cor- 
rection. 

• correct for fringes. 

• wavelength-calibrate every pixel, using the arc- 
lamp with a correction from sky-lines. 

• perform a t wo-dimensional sky-subtraction (see 
Kelsonl[2003 |. 



• extract spectra without resampling in a small spa- 
tial region around the target. 

Thus the final product of the pipeline is a large set of 
pixel fluxes for every target. Each of these pixels has an 
associated flux uncertainty, a wavelength, a wavelength 
uncertainty, and the value of the target spatial profile at 
the pixel. 
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Table 2 

The star-counts profile of NGC 2419. 
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Note. — Tfie quantities and To are tfie inner and outer radii of tlic annulus of effective area A^f f tliat lias mid- value rm, [= 0.5{ri + ro)] 
and that is composed of Ntot stars at an average radius (r) with standard deviation Ur. The surface brightness, not corrected for extinction, 
is /ivi with Poisson uncertainties 5fj,v The ACS survey does not fully cover the inner region of the cluster (see Fig. 1 of DOS), so to avoid 
radial biasses, we construct the profile from data in an angular range amin to a^nax from —180° to 180° (where West = 0° and North 
= -1-90°). The final column lists the source of the data, with 1=ACS, 2=ACS-f Suprime-Cam, 3=Suprime-Cam. 



It is then possible to measure the velocity of the target 
using a Bayesian approach. We adopt a simple model 
for a normalized stellar spectrum that consists of a con- 
tinuum broken by three Gaussian absorption profiles lo- 
cated at the rest wavelengths of the three prominent Cal- 
cium II triplet lines (8498.02 A, 8542.09 A, 8662.14 A). 
The stellar continuum is estimated using a robust algo- 
rithm that iteratively filters out absorption and emis- 
sion lines and computes a running median of the spec- 
trum. We multiply the normalized model by the con- 
tinuum function to compare to the observations. A 
Markov-Chain Monte-Carlo algorithm then finds the best 
Doppler shift and Ca II line widths given the individual 
pixel measurements, as well as the uncertainties on these 
parameters. Due to the high signal to noise of these ob- 
servations and the excellent wavelength solution of the 
arc-lamp reductions, the velocity errors determined in 
this way for our NGC 2419 targets are very low, typi- 
cally below Ikms^^. This can be achieved, despite the 
spectral resolution of ~ 42kms~^, because the S/N is 
sufficient to provide good centering constraints on the 
absorption lines. 

However, unlike fibre-fed spectrographs, slit spectro- 
graphs have the the disadvantage that small milling er- 
rors (or large astrometry errors) will cause a target to be 
offset with respect to the slit, giving rise to an artificial 
velocity offset if the effect is not corrected. Given the ge- 
ometry of DEIMOS, if a slit is erroneously milled 0"1 dis- 
placed along the dispersion direction, one would naively 
expect a 14kms~^ velocity shift; in practice though, the 
tolerance is less strict than this estimate suggests, since 
atmospheric turbulence smears out the incoming light. 
Indeed in poor seeing conditions such errors will be min- 
imized. Another way that this problem can arise is from 
misalignment of the mask. DEIMOS masks are generally 
aligned by centering bright stars (typically 4) in small 
boxes that have been milled into the mask; before start- 
ing an exposure the observer iterates on positional and 
rotational corrections until the stars are centered to the 



achievable accuracy, which is about 0"05. An error in 
the mask alignment along the dispersion direction shows 
up easily as a systematic offset in the OH sky-line wave- 
length, while a rotational error causes a radius-dependent 
shift, as would be expected. These misalignments are 
corrected for in the pipeline to high accuracy (better than 
Ikms-i). 

Nevertheless, the milling and astrometry errors are un- 
avoidable, and it is difficult to estimate their effect. One 
way to alleviate the problem is to make use of telluric 
absorption lines: the position of the lines is not depen- 
dent on the slit error, so they give a promi sing means to 
calibrate velocities (as argued by Simon fc Geha^_2007. ). 
Unfortunately, the telluric features are broad and time- 
dependent; from inspecting the effect on repeat measure- 
ments we found that the corrections that we derived were 
actually adding noise. We decided therefore not to in- 
clude a telluric absorption correction. 

The two observed masks have 12 stars in common, 
and the r.m.s. scatter between velocity differences is 
1.91 km s~^. This is encouraging, as it means that the 
individual measurements are a factor of ^ smaller. 

Due to the positioning errors, it is difficult for DEIMOS 
to provide us an accurate absolute velocity. To circum- 
vent this problem, we adopt the approach of bringing 
our measurements onto the zero-point defined by the 
High-Resolution E c helle S pectrometer (HIRES) study of 



Baumgardt et al. (|2009) hereafter B09). There are 7 
stars in common between the DEIMOS sample and the 
HIRES sample (39 stars), we find an average offset of 
Av — vhires — vdeimos = —3.06 km s^^ with an r.m.s. 
scatter of 2. 37kms~"'^. This r.m.s. scatter is produced by 
a combination of the HIRES error on each star quoted by 
B09 (Si^HiREs), the error on the DEIMOS velocities re- 
turned by our pipeline {Si^ijEiMoSi itself a combination 
of the wavelength calibration uncertainty and the pixel 
measurement uncertainties), and the milling/ astrometry 
error, which we will refer to as 6ast- We make the plau- 
sible assumption that Sast is independent of magnitude. 
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Figure 4. Parameter projections of the kinematic data set. The B09 sample is shown with open circles, while filled dots are from the 
present DEIMOS survey. The stars are marked either red, black, green, or blue according to their likely cluster membership. Blue points 
are very likely contaminants, as they lie more than 20kms~^ from the cluster mean. Red points mark those stars that are inside 2 times 
the velocity dispersion of the MOND model shown in Figure |6] and constitute the cleanest sample. Black and green points lie beyond that 
model, but inside 20kms~^ from the cluster mean, with theblack points having [Fe/H] < —1.3 and the green [Fe/H] > —1.3. Note that 
B09 do not provide metallicity information, so open circles are absent from panels 'b' and 'd'. 
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The condition: 



1 " 

-E 



jViMIRES - Vt,DEIMOS - f 
"i. HIRES ^ "i.DEIMOS ^ "ast 



= 1, 



where we are summing over the n — 7 stars in 
common between the two samples, requires that the 
miUing/astrometric error must be Sast = 2.25 km s~^. 
We add this additional uncertainty in quadrature to our 
pipeline velocity uncertainty estimates. Note however, 
that HIRES, being a slit spectrograph, also suffers from 
the same centering uncertainty. However, given its higher 
spectral resolution the corresponding velocity offset er- 
rors should be lower, and we take them to be negligible 
compared to the DEIMOS 5ast- 

In addition to the velocities, we also measure the Na I 
equivalent width (which can be a useful population dis- 
criminant), and Ca II line strengths to deriv e a metallic- 
ity in an identical way to that discussed in |Ibata et al. 
( |2005[ ). The metaUicity (on the jCarretta &: Gratton|199 7 
scale) is calculated as: [Fe/H] = -2.66 -h U.42[i;G'a 

0.64(VhB-V)], with SCa = 0.5WA8498 + l-0i;i^A8542 + 

0.6£'Wa8662, (Vhs — V) being a surface gravity correc- 
tion relative to the V-magnitude of the horizontal branch. 
We adopt the value oiYnB = 20.45 plarris 1996). Ta- 
ble [5] lists the star number, the mask number (or B09 
identification number), the position, the radial velocity 
measurement and corresponding uncertainty, the I-band 
magnitude and V — I color, the metaUicity, and the pro- 
jected radial distance R from the cluster center. Those 
stars with multiple observations are marked in the final 
column. 



2.1. Defining the kinematic sample 

The final sample comprises a total of 178 unique stars, 
of which 39 are taken from the B09 catalog. For those 
stars with multiple observations, we choose the mea- 
surement with the lowest estimated velocity uncertainty. 
Various parameter projections of the sample are shown 
in Figure |4] velocity as a function of radius, veloc- 
ity as a function of metaUicity, position on the sky, 
the metallicity-radius relation, color-magnitude position, 
and finally the velocity uncertainty-magnitude relation. 
While the majority of the observed stars are clumped 
close to the expected mean Heliocentric velocity of the 
cluster (Fig. |4^), there are clearly Galactic interlopers 
in the sample. The most obvious interlopers are marked 
blue and have velocity differences > 20kms~^ with re- 
spect to the cluster mean. In an attempt to make an 
initial selection to help guide the discussion, we take the 
velocity dispersion profile CT„(r) of the {M/L)v — 1.346 
anisotropic MOND model (#24) discussed in Section [3] 
below. This model has the highest likelihood of alT 
the MOND models we consider (for sample A), and its 
adoption in this discussion should favor the case for the 
MOND theory. Selecting stars within 2 x ct^ (r) allows us 
to make the initial separation between very likely clus- 
ter members (red) and the kinematic outliers (black and 
green) in Fig. |4^. Inspection of Figure suggests that 
a metaUicity cut at [Fe/H] = — 1.3 provides a means to 
sort the kinematic outliers into metal-poor (black) and 
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Figure 5. Spectra representative of bright, medium and faint 
stars in the spectroscopic sample that are confirmed members of 
NGC 2419. The wavelength range shown highlights the main region 
of interest around the prominent Ca II triplet lines at 8498.02, 

o 

8542.09, and 8662.14 A. The error spectrum is also shown (red 
line) in each panel, as well as the identification number of the star 
in Table [sl Note that the pipeline discussed in the text produces 
non-binned versions of the spectra (which are just pixel fluxes with 
their corresponding wavelengths); here such data have been binned 
up for display purposes only. 



metal-rich (green) samples^. Unfortunately the other 
measured parameters do not provide us with strong addi- 
tional information that could help discriminate between 
cluster members and interlopers (our Na I line strength 
cut EWNai < 2.0 rejects one dwarf star that would 
otherwise have been rejected anyway by velocity). 

Due to its importance, in Figure [6] we reproduce 
the velocity-radius relation, zooming-in on the cluster 
region, and displaying the velocity uncertainties. Note 
that we have refined the cluster mean velocity quoted by 
B09 (— 20.3±0.7kms~^), using instead a weighted mean 
of their data that gives Vmean = — 20.9kms^^; with this 
value the DEIMOS velocity distribution is noticeably 
more symmetric. The lines delineate 2 x cr(r) of the 
MOND model mentioned previously, a limit chosen to 
incorporate all of the sample of B09. Inspection of Fig- 



orp' 

ureli] shows that there are no stars with [Fe/H] < —1.3 



9 Recently, |Cohen et al.||2010| l have reported a significant metal- 
licity spread iii iNUU 2419, a property we confirm in Figure|4jl. Our 
metaUicity distribution appears wider however, probably due to the 
inclusion of lower S/N spectra from stars at the faint end of the 
color-magnitude diagram. 
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Radius [arcmin] 



Figure 6. Closeup of the cluster region of Figure |4|i, showing 
velocity uncertainties. The graph markers have the same coding as 
in that figure. The line showrs 2 X o'(r) according to the anisotropic 
{M/L)v = 1.346 MONO model (model #24 in Table|4]l. Note that 
there are no plausible cluster members outside of this it20kms~^ 
velocity window (see Figure |4^). 

that have velocities beyond 20kms^-'^ of the cluster 
mean, which in turn strongly suggests all stars with 
[Fe/H] < —1.3 are cluster members. This includes all 
stars marked in black in Figure [61 many of which have 
error bars that overlap the \v{r)\ < 2 x a{r) region. 
Nevertheless, several of the black data points lie well 
outside of that region and, at face value, greatly enhance 
the local velocity dispersion. We will return to this issue 
in Section|4] The three green points at (i?, v — Vmean) = 
(3!1, 16.8kms-i), (9^0, IS.Tkms"!), (9^0, -ll.lkms^i), 
are most probably Galactic interlopers, and can be 
safely removed from the sample. Of the four red points 
at i? > 7', it is interesting that three are metal rich, 
having -1.07 < [Fe/H] < -0.82, substantially higher 
than the cluster mean as can be seen from Fig. ^p. The 
presence of these stars is a fascinating puzzle which we 
hope to return to in a future contribution. However, 
since their velocities lie extremely close the cluster 
mean, it is inconceivable that they are not members, so 
for the purposes of the current discussion we will treat 
these stars as dynamical tracers. 

The sample whose kinematics we will analyze below 
consists of all 157 stars marked red in Figure [6| which 
we will refer to as "Sample A", and the 9 stars marked 
in black (Sample B). 



T 




J I I I 1 1 1 I I 1 1 1 I I 1 1 : 

-100 100 

PA [deg] 

Figure 7. Radial velocity as a function of position angle East 
of North for sample 'A', but removing stars in the central 0^5. 
The dashed curve shows the B09 fit to their sample. In contrast, 
fitting the data shown here gives the continuous curve, which has 
a slightly lower amplitude of 2.6kms~^, and has a 5% chance of 
being simply a spurious result due to random sampling. However, 
fitting on the DEIMOS stars alone gives the dotted fit, which is 
certainly not significant. 

2.2. Rotation 

One of the several interesting results reported by B09 
was the signature of rotation in their kinematic sam- 
ple, which they deduced had a rotational amplitude of 
3.26 ± 0.85kms-i directed towards PA = 40°9 ± 17°8, 
approximately at right angles to the direction of elonga- 
tion of the cluster found by B07. The rotation of the clus- 
ter is a critical issue, since rotational support can offset 
pressure support, and obviously influence the kinematic 
models one has to fit. In reanalyzing the B09 kinematics, 
we noticed that a significant part of the rotation signal 
was coming from a small number of stars very close to the 
cluster center. Because the dispersion increases rapidly 
towards R — 0, such stars will generally have velocities 
quite far from the cluster mean, and if the sample is small 
(as in the case of B09) stochastic effects can cause a large 
but spurious rotation signal. To examine this possibility, 
we analyze sample A, ignoring the stars with R < 0'.5; we 
find that the rotation amplitude diminishes to 2.6kms~^ 
at PA = —4° (see Figure [7|. To place confidence limits 
on the velocity amplitude we undertook a Monte-Carlo 
experiment, keeping the measured velocities and veloc- 
ity errors, while setting the position angles of the stars to 
random values. We found that 5.1% of the random tests 
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Table 3 

Velocity dispersion estimates for the two samples. 





Sample A 




Sample 


A+B 


N 


R 






R 


(T 




(arcmin) 


( km s^"*") 




(arcmin) 


(kms-i) 


26 


0.343 


6-573lJ;^i3 


27 


0.352 


'■UDi_Q ggg 


26 


0.912 


3.816tJ;0»? 


27 


0.920 


q 77Q + 1.002 
■3. 1 'O_o.443 


26 


1.340 




27 


1.355 


3.896_o 504 


26 


1.978 


9 fifi9+^ 731 


27 


2.011 


o.iuy_Q 4gg 


26 


2.613 


n 007+0.576 
u.»»l -0.549 


27 


2.625 


1 004+0.664 


27 


5.068 


75fi+0-436 


31 


5.226 


1 9Qc:+0.683 
l.Z30_o 3g-^ 



gave fake rotation amplitudes as large as that measured. 
However, if we limit ourselves to only the DEIMOS sam- 
ple, the derived rotation amplitude drops to 0.6kms~^ 
at PA — — 86°9, yet 39% of random realizations give 
rise to amplitudes greater than this value. We conclude, 
therefore, that the B09 rotation measure was largely due 
to low-number statistics, and that the evidence for sig- 
nificant rotation of the cluster is marginal. 

2.3. Velocity dispersion profile 

While the full likelihood analysis that we present in sec- 
tion His clearly the preferred and most powerful method 
for ruling out the dynamical models (and their underly- 
ing theories of gravity), it may be useful to some readers 
to visualize the velocity dispersion profile. To this end we 
calculate the estimates of the velocity dispersion, using a 
maximum-likelihood method, under the assumption that 
the underlying distribu tions are Gaussian (see Eqn. 3 of 
Pryor & Meylan 1993). For sample A, we truncate the 
Gaussian distributions to account for the effect of the 
imposed selection function. The resulting dispersions for 
samples A and A-f B are shown in Fig. [8] and reported 
in Table [3] We expect reality to lie somewhere between 
these two estimates. 

3. MODELS 

3.1. Dynamical models 

The models we compare to in this paper have been 
specificall y constructed as an ex tension of those de- 
scribed in SoUima & Nipotil pOTol hereafter SNIO). The 



phase-space distribution ot stars is given by the distribu- 
tion function 



fiE,L) = fo cxp - 



2(7 



K' a 



exp 



E 



- 1 



f{r,vr,vt) = fo exp 



exp 



2o\ 



- 1 



(2) 



( Michie [1963 1 , where E and L are, respectively, the en- 
ergy and angular momentum per unit mass, Vr and Vt are 
the radial and tangential components of the velocity, the 
effective potential V' is the difference between the clus- 
ter potential at a given radius r and the potential at the 



cluster tidal radius tp — ((> — 4>t, fo is a scale factor, and 
ax is a normalization term which is proportional to the 
central velocity dispersion. As usual, the above distribu- 
tion function has been integrated to obtain the density 
and the radial and tangential components of the velocity 
dispersion: 



2il> 1-^-2^— vl 



An 

Pir) Jo 
P{r) Jo Jo 



Vt f{r,Vr,Vt) dvtdvr 



Vt fir,Vr,Vt) dVtdVr, 



vf f{r,Vr,Vt) dvtdvr. 



The above equations can be written in terms of dimen- 
sionless quantities by substituting 



c = 



2ay 



2ai 



W 



P_ 
2 ' 



where po is the central cluster density and 
_/ ^aK \" 



is the core radius (King 1966). The parameter fa deter- 
mines the radius at which orbits become more radially 
biassed. As — >• oo, models become isotropic. The 
dimensionless expressions are: 



p{f) = 



r''^e"^(e^-''-C-l) dry dC 



S^xlCC^C-^e-f (e— C 



-^l^e^o-n-C - 1) d?7 dC 



1) dry dC 



Jo C 5 Jo e 



w-c 



1) dry dC 

1) d77 dC 



Jo C Mo 



1 fW-c 



1) dr] dC 



(3) 



The above equations allow one to derive the density and 
the velocity dispersions once the potential profile VF(f) is 
known. This last quantity is linked to the density by the 
Poisson equation (or its MOND modification; Eqn. [I]) 
with the boundary conditions at the centre 



W = Wo, 
dW 
df 



-0. 



Here (and below, when not specified otherwise) we adopt 
the "simple" in terpolating function fi{x) = x/{l + x) 
( Famaey &: Binncy 1^2005, ) . In spherical symmetry Equa- 
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Figure 8. The line of siglit velocity dispersion profile of the 26 dynamical models is overlaid on the observed profile, derived from sample 
'A' (green) and sample 'A+B' (black). Newtonian models are colored blue, MOND models red, while isotropic models are represented with 
continuous lines and anisotropic models with dashed lines. As discussed in the text, we expect reality to lie between samples A and A+B; 
the closest fits are clearly the anisotropic Newtonian models. 



tion ([T]) can be written in the form 



1 d 

j,2 f^j. 



dr 



dtp 
dr 



or, using dimcnsionless quantities, 
1 d 

f2 f-ip 





dW 


\ dW 




dr 


J 



= -9p, 



(4) 



where x = (^'k/o.qI'c is a dimcnsionless parameter, which 
is smaUer for systems closer to the deep-MOND regime 
(see SNIO). For a given choice of (Wo,Xi^o), Equations 
[3] and |4] have been integrated to obtain the 3D density 
and the radial and tangential components of the velocity 
dispersion. As a last step, the above profiles have been 
projected on the plane of the sky to obtain the surface 
mass density 



prdf 



and the line-of-sight velocity dispersion 



<7l{R) 



1 



a^R^] df 



7.^ f2 _ J^2 



(6) 



The described procedure is straightforward and pro- 
duces spherically symmetric self-consistent equilibrium 
models, truncated at a tidal radius ft. Note that spheri- 
cal symmetry is justified as the evidence to the contrary 
(a measurement of e = 0.19 ± 0.15 or e = 0.14 ± 0.07 
by B07 — depending on the sample choice) is not very 
strong. Also, the above models assume a single distri- 
bution function for all stars, regardless of their masses. 
This can be considered a good approximation since mass 
segregation effects are expected to be negligible in NGC 
2419. This can be deduced by the fact that the relax- 
ation time is s ignificantly larger than the cluster age 



3.5; Harris 1996 Marm-Franch et al. 20091 



and it is confirmed by the lack ot radial segregation of 



(5) the BSS population (DOS). 
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As will be described in Sect. |4j the full likelihood anal- 
ysis performed here requires for each model j the knowl- 
edge of the distribution of projected velocities (instead 
of simply its second-order moment). This function, here- 
after referred to as f^{v, R), can be calculated as 



r{v,R) 



p{r)r 



Vr2 - i?2 Jo 

f{r,KiVt) dvydv^dr, (7) 



adopting the following change of variables 
vl^v sin a + Vx cos a, 
v^ — [{v cosa — Va: siu Q;)^ -|- w^] , 

where 

a = arccos(i?/r) . 

3.2. Stability analysis 

It is well known that extremely radially anisotropic 
Newtonian models are pro n e to bar instability (e .g. , iFrid-i 



man fc Poliachenko 
cently shown that 



19841. Nipoti et al. pOll) have re 



his radial-orbit instability operates 



pc 



also in MOND, thus limiting the fraction of radial orbits 
in realistic models. To test the stability of our models 
we run a set of N-bod y simulations using the MOND 



A^-bod y code n-mody (Nipoti et al. 
Nipoti,^2009), which can be used to 



2007a 
oliow 



Londrillo fc 
re evolution 



of either MOND or Newtonian collisionless stellar sys- 
tems. The use of a collisionless N-body code is justi- 
fied by th e fa c t th at NGC 2419 has a long relaxation 
time (see ^3.1 [3.3 a nd BOS). The stabili ty analysis was 
performed tollowmg Nipoti et al. (2011), and we refer 
the reader to that paper tor details! The initial condi- 
tions of the simulations consist of an N-body system of 
8 X 10^ particles, which are distributed in phase space 
with a standard rejection technique usingthe analytic 
distribution function given in Equation ([2]). We veri- 
fied that the isotropic or mildly anisotropic (i.e. with 
large anisotropy radius as compared to the core ra- 
dius Tc) MOND and Newtonian N-body models are sta- 
ble. Then we explored the behavior of the models for de- 
creasing Va , finding in each case the minimum anisotropy 
radius for stability ras- For Newtonian models we found 



? = £,hai f — ^ ^OT isotropic systems) . The maximum 
values for stabihty are 1-9 and ^haif,s ^ 1-5 in New- 
tonian gravity, and ~ 2 — 2.1 and ^haif s ^ 1-4 — 1.5 
in MOND. 

3.3. Binaries 

Before proceeding to the detailed comparison of the 
above models to the new data presented in this contri- 
bution, we first consider the likely effect of binary stars 
in the sample on the distribution o f observed kinemat- 
ics. Follow ing the seminal paper by 'Rubenstein &: BaT] 



lyn (1997 hereafter RB97), the most comprehensive re- 
cent attempts to determine the binary fraction (bf) in 
globular clusters have been performed by modeling the 
spread of the Main Sequence in very accurate Hubble 
Space Telescope color-magnitude diagrams, taking into 
account the observational effects (in particular blending) 
by means of extensive art ificial star experiments ( Bcllazz^ 
ini et al.|[2002{ |Sollima et al., 2007; ,Milone et al. 2{)08f 
I'hese results are clearly dependent on the assumed dis- 
tribution of mass ratios q; however, essentially all studies 
try to overcome this problem by using several different 
distributions. 

From these analyses globular clusters have been found 
to contain copious quantities of binaries. RB97 found 
15% < 6/ < 38% in the core of NGC 6752, which is 
a much denser cluster than NGC 2419 (with logarithm 
of central luminosity density logpo = 4.9L0/pc'^ as op- 
posed to 1.5L(Tj/pc'^), and bj < 16% outside. Similarly, 



Bellazziniet al.| ( [2002| ) found 10% < 6/ < 20% in the core 
of NGC 288 (which with logpo = 1.8L0/pc^ is closer in 
density to NGC 2419) and < 10% outside the core (con- 
sistent w ith 0%). With a larg er sample of 13 low density 
clusters, Sollima et al. ( |2007 ) concluded that the 



mum bf m the central areas they sampled was > 6% in 
all the considered clusters, but may be as large as 20% in 
some cases. Most of their best fit values range between 
10% and 20%, with peaks of 50% (Terzan 7). The clus- 
ter with the central density most similar to NGC 2419 in 
their sample (NGC 6101) has a minimum bf — 9% and 
best fit values ranging from 15% to 21%, depending on 
the assumed distribution of q. 

We know that binaries must also be present in 
NGC2419, since th e cluster has numerous B lue Strag- 
glers Stars (BSS) (Dalessandro et al. 2008 hereafter 



s/vc ^ 0.9. For MOND models the value of Vas/rc can D08) while collisions are nearly impossible (because of 



depend on the mass-to- light ratio {M/L)y, so for each 
value of {M/L)v we ran a set of MOND simulations with 
decreasing Va/rc, thus determining ras/rc- The most 
anisotropic stable MOND Michie models, for the consid- 
ered {M/L)v are those reported in Tablellj with Vas/rc 
in the range 1.4 — 1.6. In terms of the half-mass radius 
Thaif, we have /r half ^ 0-5 in the Newtonian case and 
r as n half ^ 0.8 — 1 in MOND: thus Newtonian models 
of NGC 2419 can have significantly smaller ras than cor- 
respon ding MOND models , in agreement with the find- 
ings of Nipoti et al. (2011) for a different class of stellar 
systems. It is also interesting to consider the Fridman- 
Polyachenko-Shukhman parameter ^ = 2TrlTt (where 
Tr and Tt are the radial and tangential components of 
the k inetic energy tensor; see |Fridman fc Poliachenko| 
1984) and of the quantity £,haif , defined as twice the ra- 
tio or radial to tangential kinetic energy within rhai f (so 



the long relaxation time), thus essentially all BSS must 
come from the evolution of Primordial Binaries. As a 
comparison, the fraction of BSS normalized to the total 
luminosity is double that in a;-C en (whose binary fra ction 
has been estimated to be 13%, Sollima et al.|[2007 ). On 
the other hand, the ratio of BSS to Horizontal Branch 
stars Nbss/^hb = 0-3- The same par a.meter is 0.07 
in M 13, 0.3 in M3 and 0.92 in NGC 288 ( [Ferraro et aL 
2003 computed in a homogeneous way). As already re- 
ported in Sect. [3j DOS found that the BSS show no sign 
of radial segregation, supporting the idea that two-body 
relaxation is highly inefficient in NGC 2419. This is the 
observational foundation for two key assumptions of our 
analysis: first, this implies that 6/ is constant with ra- 
dius, making it easier to model a population of binaries; 
second, and more importantly, it implies that the mass to 
light ratio of the cluster is constant with radius. If mass 
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Figure 9. The star counts profile derived in jj2|(and listed in Table [2| is marked by black dots with their associated uncertainties. The 
blue dashed line on the left-hand panel shows tne best-fit model (the anisotropic Newtonian model #17). The green crosses mark the 
values of this model integrated over the same radial bin as each observed data point, taking into account the complex spatial selection 
function (gaps between CCDs, etc). The right-hand panel displays the full library of Michie models; the colors and line styles have the 
same meaning as in Fig. [s] 



segregation had been at work, as in most other clusters, a 
larger fraction of low-mass stars would have been present 
in the outer regions, enhancing [M/L)v there with re- 
s pect to the cor e. 



Sollima et al. (2008) found that the fraction F of BSS 



to mam-sequence stars (computed in a way that may be 
difficult to apply to NGC 2419) correlates very well with 
6/, and My, in the sense that brighter clusters have a 
smaller and, as a consequence, bf. So, while there 
is a possibility that bf may be lower than average in 
NGC 2419, a fraction 6/ as large as 20% is clearly possible 
and cannot be discounted. 

Note that the effect of binaries on the velocity disper- 
sion should be larger in the outer parts of the system, 
enhancing the velocity dispersion as would be expected 
if MOND were correct. 



To construct a model for the effect of the binary popu- 
lation on the observed kinematics, we run a Monte Carlo 
simulation. Following McConnachie & Cote (2010) the 
projected motion of the primary component in a binary 
system is given by 



27rai sini 

V = ^TTT „, . ,„ icos 



P(l-e2)i/2 



[cos(6' + + e cos u] 



where Oi is the semi-major axis of the primary compo- 
nent, P is the orbital period, e the eccentricity, i the in- 
clination angle to the line-of-sight, 9 the phase from the 
periastron and uj the longitude of the periastron. The 
semi-major axis has been calculated using the third Ke- 
pler law 



3.4. Modeling the binary fraction 



ai 
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P2c;(mi+m2) 
47r2 
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Figure 10. The surface brightness residuals from the model (#17) that fits best the surface brightness profile. The identification numbers 
of the models are indicated. The meaning of the dots, crosses, colors and line styles are the same as in Fig. |9] 



where mi and m2 are the masses of the primary and 
secondary component. For each simulated binary we ex- 
ctracted randomly a combination of (mi, m2, P, e, 0, u), i) 
from suitable distributions and derived the correspond- 
ing projected velocity v. We adopted a fixed mass for 
the primary component of mi — 0.83 Mq (i.e. the typ- 
ical mass of a RGB star, calculated by comparing the 
color-magnitude diagram of NGC 2419 with a suitable 
isochrone of Marigo et al. 2008) and extracted the mass 
of the secondary component from the mass-ratio distri- 
bution by Fisher ct al. (200^. We followed the prescrip- 
tions of [Puquennoy fc Mayor (1991) for the distribution 
of periods and eccentricities. We removed all those bina- 
ries whose corresponding semi-axis lie outside the range 
o-min < a < 100 AU where Omm is linked to the radius 



of th e secondary component (according to Nelson et al. 
19861. The distribution of the angles {i,6,uj) has been 



chosen according to the corresponding probability distri 
butions (Prob(i) cx sin i; Prob(6') cx 9~^; Prob(a; 
constant). 

The resulting distribution g(v) is shown in Figure 



The velocity distribution corrected for binaries fl{v, Ri) 
is derived by convolving g{v) with the dynamical mod- 
els. Since &/ is expected not to change significantly with 
radius in NGC 2419, as discussed above, we simply take: 

fl{v,R,) = {l~bf)f{v,R,)+bf [ P{v',RM^' -v)dv' , 
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(8) 

where Ri is the projected radius of star i in the sam- 
ple. Note that although the wings of the distribution in 
Figure [TT] extend to large velocities, the central peak at 
u = is very prominent, so the effect of convolving this 
distribution with the dynamical models should be small 
if the binary fraction is small. 

4. ANALYSIS 

The final step in constructing the line of sight velocity 
distribution models is to convolve f^{v,Ri) by a Gaus- 
sian function of dispersion equal to the estimated veloc- 
ity uncertainty Svi of star i. This results in the predicted 
probability distribution f^{v, Ri). The likelihood of each 
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Figure 11. The velocity distribution model for binaries. To ac- 
count for a possible binary fraction bf, this distribution is con- 
volved with the line of sight distribution of the dynamical models 
as detailed in Eqn. |8] 
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Figure 12. The anisotropy profile for model #17. 



surements at the radii of the kinematic data, what is the 
hkehhood that the data were drawn from the same distri- 
bution as the model?" . If instead we had adopted a fit 
to the velocity dispersion measures, we would effectively 
have been assuming that the velocity data sample the 
radial bins in a fair way. That is not a good assumption. 

To assess the effect of contamination in the stellar sam- 
ple, we perform the likelihood analysis on the full data set 
(samples A-l-B), but also on the more restrictive sample 
A. For this latter case, we implement a window function 
in the likelihood calculation, curtailing the wings of the 
model velocity distributions to within the selection 
region of the sample (see Figure [6]) . The model distribu- 
tions are re-normalised to compensate for the discarded 
velocity parameter space. Table|4]lists the model number 
J, the adopted theory of gravity, the interpolation func- 
tion fj, (for MOND), the assumed Heliocentric distance 
d, and the structural parameters of the cluster model. 
These parameters are the anisotropy radius To, the core 
radius r^, the half-mass radius r^, the tidal radius rj, and 
the central line of sight velocity dispersion ctq- The (base- 
10) logarithm of the ratio of likelihoods of model j to the 
best model is listed in the columns marked log ^^"-^ , both 
for samples A and A-l-B. We also list the corresponding 
{M/L)v (in Solar units); for the Newtonian case, the 
quoted values are the best-fit results, while for MOND 
they correspond to the input (fixed) value that was sim- 
ulated. The derived binary fraction &/ is also given. 

Inspection of Table |4] shows #17 to be the best rep- 
resentation of the cluster among these Michie models. 
The prefer red (M/L)v transpires to be identica l to that 
derived by McLaughlin & van der Marel (2005). To in- 
vestigate whether this best titting dynamical model is an 
acceptable representation, we undertook a Monte-Carlo 
test, drawing artificial stars from the model and com- 
puting the sample likelihood. Each of the 10000 arti- 
ficial samples is required to have the same number of 
stars as the actual kinematic sample, and the artificial 
stars are selected at the model radii of the observed stars. 
We find that 61% of these randomly-drawn samples have 
likelihood larger than that measured for model #17 from 
the observed sample. This indicates that the best-fitting 
Michie model is an excellent representation of the clus- 
ter 



10 



dynamical model j can now finally be calculated via: 

i—l 

where we have substituted the actual observed radial ve- 
locity measurements Vi into the probability distribution. 
The second sum in Equation |9] is performed over the 
surface brightness measurements E for which the model 
predictions of each model j were displayed above in 
Fig. [9] In calculating the likelihood of the surface bright- 
ness predictions we obviously are assuming a Gaussian 
model with dispersion equal to the estimated measure- 
ment uncertainties 6Y1. Note that the radial positions of 
the spectroscopic stars are not drawn in a fair manner, so 
their positional information cannot be used in the likeli- 
hood calculation. The first term on the right hand side 
of Eqn. [9] answers the question "given the velocity mea- 



5. DISCUSSION 

The analysis presented above shows clearly that the 
massive halo globular cluster NGC 2419, which is suffi- 
ciently distant not to be significantly affected by the ex- 
ternal field of the Milky Way, and is sufficiently extended 
so that a substantial fraction of its member stars expe- 
rience extremely low accelerations, conforms to our ex- 
pectations of a plausible dynamical system in Newtonian 
gravity. The Michie model #17 with {M/L)v = 1.903 
and Ta = 1.5rc, which has the anisotropy profile shown 



in Fig. 12 is found to give an excellent representation 
of the cluster. Isotropic Michie (i.e. King) models in 
Newtonian gravity are a factor of 10^^ less likely than 

Had we found a very high fraction (say 90%) of random real- 
izations with higher likelihood than that measured, it would have 
cast doubt on the model, while a very low fraction (say 10%) would 
have called the reliability of the uncertainty estimates into ques- 
tion. 
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Table 4 

The cluster structural parameters, and the adopted theory of gravity are listed for the 26 Michie models. The final six columns show the 
corresponding best fits given the data in sample A and sample A+B. The likelihoods are listed as ratios with respect to the best-fitting 
model (for each sample separately). The best-fit binary fractions are also given. Thus for sample A, model #1 is lo^^^-^ less likely than 
model #17. For the Newtonian models the {M/L)v listed are best-fit values, whereas in MOND they correspond to an input of the 

simulation. 
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the best anisotropic model, and can be firmly excluded. 
The isotropic King models in MOND for any plausible 
mass to light ratio are also completely excluded, as their 
likelihoods are exceedingly low compared to the best 
anisotropic Newtonian model. The anisotropic Michie 
models in MOND fare better: for sample A, the most 
likely such model is #26, but this is still a factor of 
40000 less likely than model #17, while for the full sam- 
ple A-l-B, the best MOND model is #24, but again it is 
much less likely (by a factor of 2 x 10'°^) than model #17. 
Therefore we also reject the anisotropic Michie models in 
MOND. 



Stellar popul atio n models by Maraston (19981, Perci- 
val et al. ( 2009 ) and Buzzoni ( 1989 ) , with age and metal- 



licity similar to NGC 2419, predict a mass to light ratio 
in the range (M/L)y = 1.6 - 3.0, dependent on the as- 
sumed IMF and to a lesser extent, on the particular evo- 
lutionary synthesis model. Values of {M/L)v ^ 3 result 
from a Salpeter IMF, while values of {M/L)y ;$ 2 arise 
from assuming a Kroupa IMF. Our best-fit model #17 
has a mass to light ratio that seems perfectly plausible 
in this context. 

Sample A was constructed to avoid velocity outliers, so 
it is natural that the MCMC algorithm returned 6/ = 
for that set. For the full sample A-l-B, the most likely 
Newtonian models have bf ^ 8%, in r easonable agree- 



ment with the findings of Sollima et al. (2007). The best 
MOND solutions however, favor bf = 0, clearly higher b / 
would increase the effective velocity dispersion, produc- 
ing a more acute discrepancy with the kinematic data. 
It is interesting to note that including a plausible bi- 
nary fraction makes the situation significantly worse for 
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Figure 13. The line of sight velocity dispersion profile of two ex- 
treme kinematic toy models of the cluster in MOND are compared 
to the dynamical Newtonian Michie model #17 (dashed blue line). 
A purely tangential model (magenta dotted line) is inconsistent 
with the data beyond the central regions, but a maximally radial 
model (full turquoise line) appears more promising. The data are 
represented with the same coding as in Fig. [s] 

MOND. 

The analysis above has considered a realistic family of 
dynamical models — the Michie models — which are con- 
sistent in the sense th at th ey have a positive distribution 
function (see Section 3.1). Michie models are generic 
models that fit real systems in Newtonian gravity (in- 
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eluding this cluster) , and refleet the expectation that the 
outer regions retain some memory of their initial col- 
lapse with higher radial than tangential motions. Alter- 
natively, one mechanism that could give rise to such ra- 
dial anisotropy, is if the cluster expan ds violently via gas 



expul sion shortly after its formation ( [Baumgardt et ciL 
20081). The considered models were also tested tor sta- 



bility, as detailed above. Although the Michie models 
are fairly general, the cluster could of course have some 
other anisotropy profile, but the difficulty is finding an 
appropriate self-consistent dynamical model that fits the 
system. 

It is perhaps instructive first to see what anisotropy 
profile would help bring the line of sight velocity disper- 
sion profile closer to the data, withou t w orrying about 
the dynamical consequences. Figure |13| compares the 
binned data previously displayed in Fig. 181 with two toy 
models. Here we have taken the density profile of the 
Michie model #17 as a given, and assume a constant 
mass to light ratio, as supported by the results of DOB. 
The magenta dotted line displays a model in MOND 
(with {M/L)v = 1.9) which has perfectly tangential or- 
bits. While this approximates the value of the central line 
of sight velocity dispersion, it vastly over-predicts that at 
larger radius. In contrast, a highly radially anisotropic 
model in MOND, that respects the Global Density Slope - 
Anisotropy Inequality (GDSAI,[Ciotti fc Morganti|2010"|, 



2/3(r) < 



din p{r) 
d ln(r) 



(10) 



which is a requirement for consistency of a wide range of 
models, is shown with the full turquoise line (for the par- 
ticular choice of {M/L)v — 0.5). Here, p is the density 

and /? = 1 — v^/v"^ is the anisotropy parameter. Clearly, 
the choice of low {M/L)y and high radial anisotropy is 
more promising. However, we recall that the GDSAI is 
a necessary (but not sufficient) condition for consistency, 
so we have no guarantee that there is a consistent model 
with such an extreme anisotropy profile, not to mention 
the issue of stability and the astrophysical inconsistency 
of such a low M/L. 

Before attempting any further modelling it will be 
helpful to assess the shape of the velocity distribution 
that the data display. The difficulty with this is that 
the amount of kinematic data is rather limited, so it is 
not very informative to show directly the observed line 
of sight velocity distribution at various radii. So instead 
in Fig. [T4j we show the predicted line of sight velocity 
distributions of the Michie model #17 at the radii of the 
kinematic observations. The distributions for the inner 
half of the kinematic data are shown in black, the outer 
half in red. These distributions have been scaled to have 
r.m.s. dispersion of unity. Note that despite the fact 
that the model has a strongly varying anisotropy (see 
Fig. 12), the integration along the line of sight yields 
distributions of very similar Gaussian shape at each ra- 
dius. The green histogram shows the distribution of ob- 
served radial velocities, where each datum is divided by 
the r.m.s. dispersion of the line of sight velocity distri- 
bution of the Michie model at the corresponding radius. 
The data clearly follow closely a Gaussian distribution 
(obviously the precise statistic depends on how one deals 
with outliers, but as we have shown above, the Michie 




Figure 14. The line of sight velocity distributions of the Michie 
model #17 that fits the cluster are shown for each radial position 
at which we have a kinematic measurement. The inner half of the 
sample is shown in black, the outer half in red. These distributions 
are shown re-scaled so that the root mean square (r.m.s.) disper- 
sion of the distributions is unity, which highlights the fact that 
they are very similar throughout the radial range of our data. The 
observed velocity distribution is also displayed (green histogram), 
with each point re-scaled by the r.m.s. of the corresponding model. 
Both the observed and best-fitting model distributions suggest that 
Gaussian line of sight velocity distributions are required. 

model #17 is an excellent fit to the data). Any model 
that we consider has to conform to this fact. 

5.1. MCMC Jeans equation analysis 

T he h ighly radially- anisotropic toy model shown in 
Fig. [Is] (turquoise line), suggests that a more general 
modcTmay help to relieve the pressure on MOND set 
by our data. The option we choose to explore next is 
to use the Jeans equation to determine profiles of the 
radial and tangential velocity dispersion that are consis- 
tent with the density profile and potential. The Jeans 



equation for a sph erically-symmetric system is (Binney 
fc Tremaine][l987| Eqn 4-55): 



GM{r) 



dlnp ^ dlnw^ ^ ^ / 
din?' dlnr \ 7H 



(11) 



where M{r) is the cumulative mass inside a radius r. 

We assume that the density profile conforms accurately 
to that of the Michie model #17, and fix the potential by 
selecting a value for the mass to light ratio. We choose 
to simulate the set {M/L)v = {1.6, 1.8, . . . , 3}, which 
span the astrophysically-plausible range according to the 
evolutionary synthesis models cited above. The challenge 
is to find now functions w^(r) and v'g{r) that are allowed 
by Eqn. [TT] and are also consistent with the observed 
data. The approach we take is to employ a Markov-Chain 
Monte Carlo (MCMC) scheme to fit trial distributi ons of 
in a similar (but not identical) fashion to 



et al. 



(2009 
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We parametrise v'^lr) at 7 radial locations, 
and fit a cubic spline to these points, with the constraints 
that is positive and falls to zero at the tidal radius. 
Starting from an initial guess for the 6 free parameters, 
the algorithm solves for Vg (r) , and uses the density profile 
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Figure 15. The density of solutions of /9(r) derived using the 
Jeans equation in a Markov-Chain Monte Carlo simulation with 
10^ steps that was implemented to fit the observed velocity 
data. The underlying density model is that of the Michie model 
#17, with a MOND potential calculated taking (M/L)v = 1.6. 
Note that most of these solutions likely do not correspond to 
dynamically-possible models. 

to project these quantities onto the hne of sight at each 
radial point for which we have kinematic data. These hne 
of sight velocity distributions are then convolved with the 
uncertainty distribution appropriate for each data point. 
A binary fraction of zero is assumed. The likelihood of 
the trial model is computed from Eqn. [9] 

The MCMC algorithm we developed to sample the 
posterior probability distribution function is a new hy- 
brid scheme which we intend to detail in a forthcoming 
contribution. To sample the full para meter space th e 
algorithm employs Parallel Tempering d Gregory 2005), 
a technique similar to simulated tempering, which uses 
parallel chains to probe the likelihood surface at differ- 
ent "temperatures" (a higher temperature produces a 
shallower likelihood surface, allowing the algorithm to 
reach distant regions beyond a local maximum). How- 
ever, instead of stepping though parameter space with 
the traditional Metropolis-Hastings c riterion, we use the 
Affine-Invari ant Ensemble S ampler o f|Goodman fc Weare| 



p]10). The [Goodman fc W care (20To|) algorithm with" 
Wretch" moves is very convenient as TTnaturally adapts 
the proposal step size (and, although not of paramount 
importance to the present problem, the afhne-invariance 
property ensures that one does not need to know the rel- 
ative scales of the parameters in advance). Tests show 
that our new hybrid scheme is very powerful, outperform 



(20101 



engmg 



ing Parallel Tempering or the Goodman & Weare 
algorithm on a variety of computationally-chal 
problems. 

The particular implementation we constructed has 6 
parallel MCMC chains, each sampled with an ensemble 
of 100 points per iteration. The chains are swapped (i.e. 
"temperered" ) once every 100 iterations on average. We 
add two astrophysical priors to the solution: first, the 
solution must conform to the GDSAI (Eqn. 10), and sec- 
ond, we require that $,haif,s < 1-5, which as aiscussed in 
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Figure 16. The velocity dispersion profiles of the most likely 
model i n M OND resulting from the Jeans equation solutions shown 
in Fig. |15| The diagram demonstrates the working of the fitting 
procedure: the large dots mark the 7 points (the two inner points 
are close and appear merged) that were used to anchor the spline 
fit (dot-dashed line) to (Tr(r). Assuming the fit to the luminosity 
profile made for model #17, and with {M/L)v = 1-6, the Jeans 
equation is used to recover the crg{r) profile (dashed line). Pro- 
jecting these along the line of sight results in ctv{R) (continuous 
line). The crosses mark (t„ at the positions of the kinematic mea- 
surements after convolving with the corresponding observational 
uncertainties. A visual comparison to Fig. [6] shows clearly that 
even this model, which is the most likely of tne MOND solutions, 
over-predicts the velocity dispersion at large radius. 



§ |3.2[ avoids the most extremely radially unstable con- 
figurations. We further require that v^ir) > 0, which 
for practical (technical) reasons we implement as a very 
strong prior. Since we have found that the observed ve- 
locity distribution is closely Gaussian (Fig. 14), we take 

v^{r) and Vg{r) to be truncated Gaussians. The trun- 



cation velocity at r is set t o \/2\l/(r), where \['(r) is t he 



local relative potential (see Binney fc Tremaine||1987 ). 

Figure 15 shows 10^ steps of the Markov Chain run- 
ning in a MOND potential with {M /L)v — 1-6, the mass 
to light value that gives the highest likelihood of all the 
simulations. The /3 parameter is calculated as a func- 
tion of (non-projected) radius r from spline fits to the 
radial and tangential velocity dispersion profiles. The 
density of solutions in the Markov Chain is proportional 
to their likelihood, so it is clear from this figure that the 
data — even in MOND — favor an anisotropy profile 
that is similar to a Michie model in the central regions, 
being isotropic in the very center and mostly radially 
anisotropic out to w 5'. Between ss 5' and « 11' the so- 
lutions are highly degenerate and the data do not appear 
to favor any particular model. We have checked that 10^ 
iterations are sufficient for convergence: repeated MCMC 
simulations with very different initial parameters recover 
identical best-fit parameters to better than 0.3%. We 

are aware that our choice of anchoring Wg(r) at 7 radial 
locations (rather than some other number, and at some 
other locations), as well as the choice of employing cu- 
bic spline interpolation, and the implementation of the 
priors, likely have a significant impact on t he p robability 
density function of solutions found in Fig. [TTJ However, 
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Figure 17. External acceleration exerted by the Milky Way on 
NGC 2419 after a time interval t = 1, 5 and 10 tdyn(r) as a func- 
tion of the mass fraction comprised within r. The case of the two 
different eccentric orbits e = 0.5 and e = 1 arc shown with red and 
black lines, respectively. 

we will not explore this complex issue further, as our aim 
here is only to show that possible alternative solutions 
may exist. 

We stress that the Jeans equation analysis does not 
guarantee that the solutions are physically plausible, in- 
deed most of the anisotropy profiles in Fig. [15] probably 
do not correspond to self-consistent models (Te., models 
with positive distribution function, which are a steady- 
state solution of the coUisionless Bolzmann equation). 
Further work is required to check the solutions; and 
clearly given the number of possibilities this is a mon- 
umental task! With this caveat in mind, we note that 
the kinematics of the most likely MOND solution from 
the MCMC analysis, whose velocity dispersion profiles 
are displayed in Fig. [l6j is a factor of 350 less likely than 
the most likely solution in an identical MCMC simula- 
tion undertaken in Newtonian gravity. Although further 
work is required to rule out MOND categorically with 
the kinematics of this cluster, the present data already 
challenge the theory severely. 

5.2. External field effect 

Throughout the paper we have treated the globu- 
lar cluster NGC 2419 as isolated. This is justified by 
the large distance of the cluster from the Galactic cen- 
ter, suggesting that tidal effects are small and that, in 
MOND, the EFE is not important. However, given the 
poor performance of MOND at reproducing the kinemat- 
ics of this object, it is worth trying to estimate quanti- 
tatively the EFE and to verify whether the discrepancy 
between data and models can be accounted for by the 
EFE. As pointed out in the Introduction, the gravita- 
tional field of the Galaxy at the location of NGC 2419 
is estimated to be, in modulus, Qext ~ O.lao and di- 



rected approximately towards the Galactic center^-'^. In 
the frame of reference of the cluster, this external field 
will change in time, both in direction and in magnitude, 
as a consequence of the orbital motion of the cluster. 
The effect of such a variation on the internal kinematics 
of the cluster is linked to the ratio between the internal 
dynamical time (tdyn) and the orbital period [Porb)- Al- 
though the orbit of NGC 2419 is unknown, we can set 
an upper limit to this ratio by assuming the cluster is 
presently at its apocenter. As a first order approxima- 
tion, we calculated the orbital period assuming planar 
orbits within a spherical isothermal Galactic potential 

= 187 km s""'^ (to have 
94.8 kpc). The dynamical time 



iiMW = v^rc^i^RGC with Vc 



Qext = 0.1 oo at Rgc 
is calculated as 



idyn{r) = 27r 



9int{r) ' 



and we adopt the values of .gmt(r) predicted by the 
best MOND model #24. Following the above consid- 
erations, even assuming a free falling orbit (c=l) we 
obtain tdyn/ Porb = 0.006 at the half-mass radius and 
tdyn/ Porb = 0.123 at the tidal radius. To understand the 
effect of such a variation on the internal cluster kinemat- 
ics, we plot in Fig. [T7]the external acceleration exerted 
on the cluster after a time t as a function of the fraction of 
mass comprised within the radius where t/tdyn{f) =1, 5 
and 10, for two different eccentric orbits (e = 0.5 and 1). 
It is evident that in all cases more than 90% of the clus- 
ter feels an almost constant external acceleration (always 
smaller than 0.2 ag) after several dynamical times. So it 
is reasonable to assume that the cluster is in equilibrium 
in the MOND potential calculated with the present-day 
external field — Vtpe^t (assumed uniform). Such an equi- 
librium configuration will not be spherically symmetric 
and must be a solution of Equation ([T]) with boundary 

conditions V(/3 — )■ Vtpext for |r| — oo and density dis- 
tribution p such that, in projection, it reproduces the 
observed surface-brightness profile of NGC 2419. 

We try to build such equilibrium models using the 
MOND N-body code n-mody (see Section 3.2 ), in which 
we implemented the EFE by imposing the above bound- 
ary conditions. In practice, given a density distribu- 
tion p, the MOND field is calculated by starting from 
a guessed potential with the correct asymptotic behavior 
(V(p — ?> Vt/Jeat) and then the solution V<y9 is found with 
the same relaxatio n procedure used in the ab sence of the 
external field (see Londrillo & Nipoti 2009[ ). Of course, 
the particles are then evolved in phase-space using only 
the "internal part" of the acceleration — (V(p — Vipext)- 
While in the absence of the external field we were able 
to set up equilibrium initial N-body models using the 
distribution function ([2]), in the presence of the external 

Independent of the gravity law, a cluster located in the plane 
of the Galaxy at a distance R from the Galactic center experiences 
a gravitational field <x v'^^^{R)/R directed towards the Galactic 
centre (where Vrot is the rotational velocity of the Galaxy) . If the 
cluster is out of the plane, the field it experiences could be different 
in MOND and in Newtonian gravity with dar k matter, due to th e 
component of the field orthogonal to the plane ([Nipoti et al.|2007b| . 
However, this effect should not be strong for iNLiU ^4iy, given its 
large Galactocentric distance and relatively small distance from the 
Galactic plane. 
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Figure 18. Surface-density profiles of MOND models of 
NGC 2419 in the presence of the external field of the Galaxy and, 
for comparison, of the best-fitting Newtonian model #17. 
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Figure 19. Velocity-dispersion profiles of MOND models of 
NGC 2419 in the presence of the external field of the Galaxy and, 
for comparison, of the best-fitting Newtonian model #17. 

field we cannot do better than initializing our N-body 
simulations with quasi-cquilibrium distributions. In par- 
ticular, exploiting the fact that g^xt ^ ooj we simply 
use N-body realizations (with 8 x 10^ particles) of the 
spherically symmetric models of NGC 2419 reported in 
Table HI and let them evolve in the presence of the exter- 
nal field. We find that the models rapidly reach a new 
equilibrium configuration, which is axisymmetric with 



we plot, respectively, the circular- 
ized surface-density profiles and circularized line-of-sight 
velocity-dispersion profiles of the end-products of sim- 
ulations with g^xt = O.lao, in which the initial condi- 
tions correspond to MOND models #4 (the best-fitting 
among the isotropic models, with {M/L)v ~ 0.1), #7 
(an isotropic model with more realistic (MjV)y — 1.346) 
and #26 (the radially anisotropic model which has the 
highest likelihood among all the considered MOND mod- 
els fit to sample A). For each model we plot the profiles 
with line-of-sight both parallel and orthogonal to the ex- 
ternal field; for comparison we also plot the profiles of 
the same MOND models in the absence of the external 
field as well as that of the best-fitting Newtonian model 
#17. We note that the end-products are not highly flat- 
tened: in particular, the ellipticity is in all cases e ^ 0.01 
when the cluster is viewed along the direction of the ex- 
ternal field (which should be approximately the case for 
our observations of NGC 2419) and e ;S 0.15 when the 
line-of-sight is orthogonal to g^xt- From Fig. 18 it is 
clear that the EFE has little infiucnce on the surface 
brightness profiles: only the very low-mass model #4 
presents a significant deviation from the reference model 
#17, but only at i? ^ 7 arcmin. Similarly, the velocity 
dispersion profiles of the end-products of the EFE simula- 
tions do not deviate substantially from the pr ofile s of the 
same models with no external field (see Fig. [l9| ; again, 
the model most affected by the EFE is the {mJl)v=^.\ 
model #4, with a significantly lower velocity dispersion 
even at i? ^ 5 arcmin, where we have most spectroscopic 
data. Overall, a comparison of the MOND models with 
EFE with the best-fitting Newtonian model #17 (blue 
curve in Fig. fl9]) clearly shows that an external field of 
magnitude O.lao does not help reconcile MOND models 
with the data. We ran similar simulations with an ex- 
ternal field 5e£ct = 0.2ao, finding slightly stronger effects, 
but the same qualitative behavior. We conclude that for 
NGC 2419 the EFE is unlikely to represent the way out 
for MOND. 

6. CONCLUSIONS 

We present new spectroscopic data of the distant glob- 
ular cluster NGC 2419, the most distant of the large 
Galactic globular clusters. Taken with the DEIMOS 
instrument on the Keck II telescope, these have been 
reduced with a powerful new spectral extraction pack- 
age that was designed partly to take advantage of these 
state-of-the-art data of NGC 2419. By avoiding pixel 
resampling we are able to maintain optimal accuracy, 
and we avoid having correlated noise. The algorithm has 
yielded an accurate sample of cluster kinematics with re- 
liable velocity uncertainties (derived by an MCMC fit to 
the spectra), which have been tested against the exter- 
nal dataset presented by B09. The spectroscopic sample 
contains 166 likely cluster member stars. We have also 
reanalyzed deep and wide-field photometric star-counts 
surveys of the cluster, taken with HST/ACS and Sub- 
aru/SuprimeCam. 

A set of 26 Newtonian and MOND Michie models of 
the cluster were constructed to compare against the ob- 
servations, and N-body realizations of these models were 
evolved to ensure dynamical stability. The Michie mod- 
els were projected into velocity distributions along the 
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line of sight, convolved with the velocity uncertainties 
and a plausible model for binaries, allowing us to calcu- 
late the likelihood of each kinematic data point for each 
model. We also fold in the likelihood of the star-counts 
profile, carefully taking into account the footprint of the 
imaging surveys. Due to the easy scalability in the New- 
tonian case, it was readily feasible to scan a large range 
of {M/L)v values. For MOND, we started from plausible 
initial values and homed-in manually on the best-fitting 
{M/L)v. 

We find that isotropic Michie models in both theories 
of gravity are very much at odds with the observations, 
and can be ruled out. However, our anisotropic Newto- 
nian model with ra/vc = 1.5 and {M/L)v ~ 1.903, and 
with total mass M = 9.12 x 10^ Mq is in good statistical 
agreement with the present observations, and is there- 
fore an excellent representation of the system. We note, 
parenthetically, that this implies that the structure and 
dynamics of the cluster can be reproduced satisfactorily 
without invoking dark matter, in agreement with B09. In 
contrast, the best anisotropic MOND model is a factor 
of ~ 40000 less likely. 

To extend these conclusions to more general models, 
we implemented a Markov-Chain Monte-Carlo sampling 
routine to probe solutions to the Jeans Equation that 
are consistent with the luminosity profile and kinemat- 
ics. While such a procedure does not ensure physically- 
consistent solutions (with a non-negative distribution 
function), it nevertheless suggests that any plausible 
MOND model will have severe difficulties reproducing 
the data, as the best of the MOND solutions is still a 
factor of 350 less likely than the best Newtonian solu- 
tion. 

Finally, we show with N-body simulations that the ex- 
ternal acceleration field is unlikely to affect significantly 
the above conclusions. 

Until now, most of the existing analyses in the low 
external field regime focussed on the global veloc i ty dis - 
persion in sparse clusters (e.g. Haghi et al. 2009 2011). 
In contrast, here we make use ot the full velocity distri- 
bution and its radial dependence, thus making our test 
much less dependent on the adopted {M/L)v- Further- 
more, those authors who have studied cluster radial ve- 
locities, generally bin their kinematic data to measure 
first and second Gaussian moments. The likelihood anal- 
ysis presented here works directly on the individual ve- 
locity measurements for each star, thus conserving the 
full spatial resolution, and the full discriminative power 
of the kinematic data. 

Our conclusions rest on a number of assumptions: 



• The cluster is taken to be a spherical body; while 
we have shown that the kinematic sample does not 
display significant evidence for a preferred axis of 
rotation, there is marginal evidence for oblateness 
in the integrated surface brightness distribution at 
about the 2(T confidence level in the inner 2' (B07) . 
At larger radii the cluster appears rounder. If the 
cluster in reality departs significantly from spheri- 
cal symmetry, the present analysis will not be valid. 
Orbits within a triaxial, or indeed a more complex 
structure, will in general be very different to the 
spherical case that was modelled. Since we observe 
the cluster in projection, and have only limited 



constraints on the line of sight depth, it is pos- 
sible that the cluster is elongated or compressed 
along the line of sight. If the cluster is undergo- 
ing tidal disruption, the natural expectation is that 
the structure should be elongated in this direction 
(given that satellite orbits are predominantly ra- 
dial). The precise kinematic behavior is hard to 
predict. If a substantial number of unbound extra- 
tidal stars are present in the sample, the co nfigura- 
tion \y ould be similar to the suggestion by |Kroupa| 
(1997) regarding dwarf satellite galaxies, and this 
would have the tendency to aument the apparent 
velocity dispersion, making the situation worse for 
MOND. In contrast, if the bound cluster stars have 
been stretched along the line of sight, conservation 
of phase space density suggests that the line of sight 
velocity dispersion should drop with respect to the 
spherical case. This might provide a get-out clause 
for MOND. While it is currently impossible to de- 
termine the three-dimensional shape of NGC 2419, 
globular clusters in nearby galaxie s such as M31 
'e.g. in the outer halo sample of 'Mackey et al 



2010) do not, to our knowledge, display significant 



morphological distortions in the direction pointing 
towards their host galaxy^^. 

We assume the cluster to be an isolated system, 
with a single self-consistent dynamical tracer pop- 
ulation. At first sight this appears a poor assump- 
tion, since NGC 2419 has been proposed at vari- 
ous tim es to be the remnant of a once nucleated 
galaxy (Mackey & van den Bergh|[2005l ), or a clus 
ter ejected from the Hagittarms strearn ( Newberg 



et al. 2003), or even 



Uinescu et al. 



the Virgo stream "fUSelTi 
While these claims remsun 



|2009D . 

somewhat conjectural at present, what concerns us 
most for the present discussion is whether the pre- 
vious evolutionary history of this structure has pos- 
sibly implanted dynamical contaminants that may 
invalidate our analysis. Clearly NGC 2419 does not 
possess the complex dynamical mixture observed 
i n M54 at the heart of the Sagittarius dw arf galaxy 
(|Bellazzini et al.|2008j [Ibata et al |2009[ ) , so the ef- 
fect is not overwhelming, but is likely to be subtle. 
Furthermore, from the analysis by B09 and from 
the discussion above, it appears that there is prob- 
ably no significant amount of dark matter bound to 
the cluster (although presumably a cored mini-halo 
with a large core radius could be accommodated). 
Thus although it is difficult to rule out the possibil- 
ity that our analysis is affected by the complexities 
of the real system with its (possible) stream and 
its (possible) dark matter halo, we estimate that 
these issues are not signi ficant. On the oth er hand, 
there is ample evidence ( Cohen et al.|2010 see also 
Fig. H above) for a small ('^ 0.2 dex) S]3read in 
metalhcity, indicating that we are not dealing with 
a single stellar population. Unfortunately, the kine- 



As we have shown in fJT] in the Milky Way the only remote 
cluster that is also massive (and has many stars to trace the velocity 
profile) is NGC2419. This is not the case in M31 where there are 
several clusters of this class, and in the ELT era it will be feasible 
to repeat the present analysis for many of them. If similar results 
are found, it will be difficult to invoke special alignment conditions. 
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matic data are not sufficiently numerous at present 
to attempt a more refined analysis, for instance, 
repeating the calculations presented above for sev- 
eral sub-samples split as a function of metallicity. 
However, inspection of Fig.HJi suggests that there 
is little radial variation of the metallicity distribu- 
tion function with radius, so it would be very sur- 
prising if such an analysis changed significantly the 
present conclusions. 

• We assume that NGC 2419 is a static structure. 
This need not be the case if the cluster moves on 
a highly radial orbit. Indeed if it was once associ- 
ated to the Sagittarius dwarf or the progenitor of 
the Virgo Stellar Stream, it is plausible that it has 
a similar pericenter to those structures (~ lOkpc). 
In that case, the cluster will have experienced a pe- 
riodic, highly varying, external acceleration field. 
Although the consequences of this possibility can- 
not be fully fathomed without undertaking detailed 
N-body simulations in a realistic Galactic potential 
with a realistic orbit, we note that the dynami- 
cal time estimated above suggests that the cluster 
settles down into equilibrium on a short timescale 
compared to its orbital period. Thus for the present 
purposes of testing gravity, NGC 2419 should be ef- 
fectively in dynamical equilibrium. 

• Implicit in our work is also the assumption of a 
constant M/L throughout the cluster. Dark rem- 
nants such as black holes or neutron stars are likely 
to be unimportant in this context, since they con- 
tribute only 3% and 0.8% to the total mass bu dget, 
respectively, acco rding to the BASTI models (Per 
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cival et aL] 2009 ). However, faint main-sequence 
stars may be more problematic. If NGC 2419 ex- 
perienced strong primordial mass-segregation, such 
as is inferred in nearby young starburst clusters 
( [Stolte et al]|2005[ |2006[ ), this would influence the 
line of sight velocity dispersion profile, lowering the 
central value and producing a shallower slope. We 
are currently attempting to measure the mass func- 
tion in NGC2419 using HST star counts with the 
Wide Field and Planetary Camera 3 to address this 
problem directly (Dalessandro et al., in prepara- 
tion); however, as discussed above, the similarity 
of the Blue Straggler Star profile to that of other 
populations already appears to exclude significant 
mass-segregation in this cluster (DOS). 



The present observations and analysis provide a strong 
test of the nature of gravity in the low acceleration 
regime. With the above assumptions and caveats, we 
conclude that the Newtonian approximation to General 
Relativity provides — by far — a better representa- 
tion for the structure and dynamics of NGC 2419 than 
MOND. 
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Table 5 

Spectroscopic measurements, with the corresponding photometry. 



star mask/ID 



IFc/HJ 

-1.862 
-3.412 
-3.451 
-1.848 
-1.828 
-1.772 
-1.862 
-1.455 
-3.447 
-3.369 
-1.360 
-1.550 
-3.501 
-2.088 
-3.484 
-1.890 
-1.916 
-3.517 
-3.483 
-1.473 
-1.436 
-3.389 
-3.434 
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RA 
(J2000) 
7 38 08.90 
7 38 08.64 
7 38 08.01 
7 38 08.00 
7 38 08.29 
7 38 07.86 
7 38 08.60 
7 38 07.66 
7 38 08.80 
7 38 07.67 
7 38 07.19 
7 38 09.44 
7 38 10.07 
7 38 09.95 
7 38 09.69 
7 38 09.27 
7 38 09.27 
7 38 07.71 
7 38 10.18 
7 38 09.47 
7 38 10.42 
7 38 06.32 
7 38 05.90 
7 38 05.90 
7 38 11.30 
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7 38 10.30 
7 38 06.42 
7 38 11.58 
7 38 05.80 
7 38 07.60 
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7 38 06.89 
7 38 12.20 
7 38 04.83 
7 38 11.00 
7 38 09.57 
7 38 05.55 
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7 38 04.49 
7 38 08.75 
7 38 13.06 
7 38 09.90 
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7 38 11.98 
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7 38 06.56 
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7 38 12.77 
7 38 07.52 
7 38 12.31 
7 38 05.94 
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7 38 13.87 
7 38 11.63 
7 38 04.16 
7 38 11.15 
7 38 14.21 
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7 38 12.41 
7 38 06.42 
7 38 03.05 
7 38 03.75 
7 38 02.87 
7 38 04.90 
7 38 14.92 
7 38 03.37 
7 38 02.74 
7 38 01.69 
7 38 01.92 
7 38 01.84 
7 38 03.29 



Dec 
(.12000) 
38 52 53.1 
38 53 00.0 
38 52 53.1 
38 52 53.2 
38 53 01.5 
38 52 56.9 
38 53 04.3 
38 53 02.0 
38 52 42.9 
38 52 45.5 
38 52 57.5 
38 53 06.3 
38 52 54.3 
38 53 04.3 
38 52 40.9 
38 52 35.7 
38 52 35.7 
38 53 15.7 
38 53 09.6 
38 52 32.0 
38 52 35.9 
38 52 38.7 
38 52 52.0 
38 52 52.1 
38 53 07.1 
38 53 02.6 
38 52 25.6 
38 53 21.5 
38 53 05.8 
38 53 17.3 
38 53 34.1 
38 53 34.2 
38 53 33.5 
38 53 33.6 
38 52 51.3 
38 53 04.1 
38 53 29.7 
38 53 38.7 
38 52 23.1 
38 52 19.7 
38 53 28.3 
38 53 37.5 
38 53 18.7 
38 53 48.8 
38 53 11.1 
38 52 00.7 
38 52 11.6 
38 53 40.7 
38 52 11.0 
38 51 59.1 
38 51 59.1 
38 51 59.8 
38 52 16.6 
38 51 52.9 
38 52 09.7 
38 53 50.8 
38 53 24.9 
38 53 24.9 
38 53 36.2 
38 52 34.2 
38 51 58.8 
38 52 10.7 
38 51 53.4 
38 52 36.0 
38 52 36.1 
38 52 02.6 
38 51 49.7 
38 53 24.4 
38 53 39.6 
38 52 24.1 
38 51 50.6 
38 52 36.5 
38 52 06.0 
38 52 14.4 
38 52 54.2 
38 53 15.4 
38 52 31.1 
38 51 57.2 



( km s ) 
-21.506 
-23.690 
-10.120 

-9.294 
-12.013 
-28.078 
-19.134 
-17.968 
-33.980 
-10.220 
-22.569 
-11.755 
-16.410 
-14.820 
-16.790 
-27.394 
-29.969 
-24.270 
-15.950 
-14.065 
-22.231 
-29.270 
-27.170 
-24.912 
-20.373 

-5.148 
-30.009 
-19.325 
-12.220 
-24.109 
-23.660 
-19.461 
-16.240 
-18.856 
-26.590 
-31.310 
-17.014 
-22.810 
-20.172 
-24.320 
-29.221 
-22.160 
-20.409 
-15.343 
-24.377 
-15.560 
-20.403 
-27.474 
-25.504 
-25.409 
-24.024 
-21.484 
-23.786 
-21.132 
-26.248 
-23.490 
-19.591 
-19.566 
-14.816 
-20.229 
-17.571 
-19.954 
-18.020 
-19.990 
-22.037 
-18.242 
-28.054 
-19.784 
-19.113 
-23.630 
-20.960 
-37.390 
-23.196 
-23.532 
-16.781 
-26.520 
-14.943 
-20.900 



5v 

( km s^-'-) 
OTT 
0.430 
0.810 
2.264 
2.270 
2.365 
2.269 
3.140 
0.640 
0.430 
2.282 
2.267 
0.670 
2.319 
0.580 
2.311 
2.289 
0.660 
0.560 
2.277 
2.482 
0.440 
0.530 
2.278 
2.285 
2.272 
2.358 
2.278 
0.610 
2.282 
0.550 
2.264 
0.710 
2.274 
2.383 
2.309 
2.282 
0.580 
2.366 
2.267 
2.270 
0.630 
2.589 
2.430 
3.169 
0.530 
2.296 
2.274 
2.278 
2.282 
2.263 
2.272 
2.480 
2.281 
2.282 
0.830 
2.301 
2.274 
2.265 
2.308 
2.263 
2.290 
1.080 
0.460 
2.360 
2.273 
2.492 
2.661 
4.424 
0.400 
1.280 
2.279 
2.582 
2.285 
2.343 
0.610 
2.301 
1.290 



I 

(mag) 
17.404 
16.315 
16.104 
16.104 
16.941 
18.556 
16.837 
16.980 
16.124 
16.487 
17.362 
17.244 
15.807 
18.308 
15.924 
18.237 
18.237 
15.778 
15.923 
17.685 
19.701 
16.413 
16.230 
16.230 
17.798 
17.370 
18.918 
17.542 
15.940 
17.758 
16.082 
16.082 
15.939 
15.939 
19.287 
17.711 
17.135 
16.113 
18.746 
16.887 
16.588 
16.037 
19.865 
16.715 
17.468 
15.635 
17.992 
17.167 
16.291 
16.895 
16.895 
17.660 
19.785 
17.721 
17.506 
16.673 
16.892 
16.892 
16.121 
18.633 
15.673 
17.628 
16.774 
16.054 
16.054 
17.702 
19.735 
19.872 
17.557 
16.215 
16.937 
17.607 
19.713 
17.085 
18.806 
15.783 
17.318 
17.074 



V-I 
(mag) 
1.176 
1.339 
1.403 
1.403 
1.239 
1.056 
1.269 
1.224 
1.398 
1.327 
1.163 
1.203 
1.514 
0.960 
1.461 
1.020 
1.020 
1.484 
1.466 
1.139 
0.853 
1.324 
1.339 
1.339 
1.105 
1.146 
0.968 
1.152 
1.440 
1.123 
1.419 
1.419 
1.455 
1.455 
0.965 
1.172 
1.211 
1.421 
0.988 
1.259 
1.301 
1.432 
0.903 
1.348 
1.195 
1.645 
1.066 
1.205 
1.370 
1.256 
1.256 
1.163 
0.966 
1.064 
1.164 
1.295 
1.239 
1.239 
1.409 
1.049 
1.584 
1.165 
1.283 
1.449 
1.449 
1.133 
0.940 
0.940 
1.163 
1.375 
1.235 
1.122 
0.944 
1.234 
0.979 
1.492 
1.160 
1.233 



R 

(arcmin) 

o:o82 

0.089 
0.102 
0.103 
0.118 
0.131 
0.158 
0.203 
0.208 
0.226 
0.261 
0.262 
0.304 
0.321 
0.327 
0.353 
0.353 
0.380 
0.407 
0.425 
0.488 
0.505 
0.510 
0.510 
0.580 
0.592 
0.600 
0.602 
0.624 
0.646 
0.677 
0.679 
0.716 
0.718 
0.721 
0.732 
0.756 
0.759 
0.783 
0.792 
0.857 
0.862 
0.877 
0.900 
0.926 
0.943 
0.977 
0.986 
0.996 
1.004 
1.004 
1.019 
1.046 
1.051 
1.056 
1.057 
1.058 
1.058 
1.076 
1.099 
1.115 
1.122 
1.147 
1.153 
1.153 
1.156 
1.160 
1.171 
1.189 
1.212 
1.281 
1.285 
1.290 
1.310 
1.327 
1.327 
1.357 
1.399 



Multiple 
observations? 



Note. — For those stars that have multiple measurements, we adopt the measurement with the lowest uncertainty. 



Ibata et al. 



Table 1 
...continued... 



71 


mask 1 


7 


38 


14. 


,59 


38 


53 


41. 


,9 


-24.363 


2 


.301 


17 


.694 


1, 


.139 


-1, 


.649 


1, 


.419 


n 


72 


B09 35 


7 


38 


03, 


,55 


38 


51 


51, 


,9 


-18.450 





.980 


17 


.090 


1, 


.218 


-3 


.236 


1, 


.426 


y 


72 


mask 2 


7 


38 


03, 


,54 


38 


51 


52, 


,0 


-18.874 


2 


.269 


17 


.090 


1, 


.218 


-1, 


.386 


1, 


.426 


y 


73 


mask 1 


7 


38 


01, 


,35 


38 


53 


27, 


.7 


-22.030 


2 


.369 


18, 


.710 


1, 


.048 


-1 


.801 


1, 


.497 


n 


74 


mask 2 


7 


38 


00, 


,77 


38 


52 


58, 


,1 


-23.904 


2 


.280 


17 


.388 


1, 


.200 


-1 


.851 


1 


.507 


n 


75 


mask 2 


7 


38 


11, 


,61 


38 


51 


31, 


,6 


-21.251 


2 


.291 


18, 


.378 


1, 


.078 


-1 


.713 


1 


.514 


n 


76 


mask 1 


7 


38 


15, 


,89 


38 


53 


25, 


,4 


-21.291 


2 


.363 


18, 


.980 


1, 


.022 


-1 


.781 


1 


.523 


n 


77 


mask 1 


7 


38 


16, 


,42 


38 


53 


01, 


,0 


-12.572 


2 


.487 


19, 


.408 


0, 


.974 


-1 


.571 


1 


.543 


n 


78 


B09 7 


7 


38 


11, 


,93 


38 


54 


19, 


,4 


-22.400 


1 


.260 


17, 


.877 


1, 


.112 


-3 


.053 


1 


.558 


n 


79 


mask 1 


7 


38 


09, 


,61 


38 


51 


21, 


,8 


-23.079 


2 


.430 


19, 


.256 


1, 


.011 


-1 


.690 


1 


.566 


n 


80 


B09 1 


7 


38 


02, 


,83 


38 


54 


01, 


,7 


-21.630 


1 


.310 


17, 


.461 


1, 


.145 


-3 


.156 


1 


.569 


n 


81 


mask 1 


7 


38 


10, 


,67 


38 


51 


23, 


,5 


-22.992 


2 


.463 


19, 


.330 


1, 


.001 


-1 


.487 


1 


.580 


n 


82 


mask 2 


7 


38 


16, 


.70 


38 


52 


36, 


.1 


-21.921 


2 


.262 


16, 


.886 


1, 


.253 


-1 


.406 


1, 


.624 


n 


83 


mask 2 


7 


38 


14, 


,96 


38 


51 


51, 


.6 


-18.102 


2 


.259 


15, 


.951 


1, 


.455 


-1, 


.856 


1, 


.640 


n 


84 


mask 1 


7 


38 


13, 


,39 


38 


54 


19, 


.2 


-26.983 


2 


.586 


19, 


.627 


0, 


.985 


-1, 


.439 


1, 


.696 


n 


85 


mask 2 


7 


37 


59, 


,90 


38 


53 


13, 


.2 


-18.270 


2 


.394 


19, 


.105 


0, 


.993 


-1, 


.634 


1, 


.703 


n 


86 


B09 SIO 


7 


38 


16, 


.92 


38 


53 


35, 


.0 


-20.750 





.580 


16, 


.227 


1, 


.384 


-3 


.423 


1, 


.768 


n 


87 


mask 2 


7 


37 


59, 


,14 


38 


53 


09, 


,3 


-17.711 


2 


.358 


18 


.715 


1, 


.044 


-1 


.483 


1 


.839 


n 


88 


mask 1 


7 


38 


02, 


,33 


38 


51 


31, 


,0 


-24.023 


2 


.275 


17, 


.292 


1, 


.201 


-1 


.966 


1 


.845 


n 


89 


mask 1 


7 


38 


17 


,46 


38 


52 


04, 


,5 


-18.671 


2 


.266 


17, 


.229 


1, 


.180 


-2 


.035 


1 


.934 


n 


90 


mask 2 


7 


38 


15, 


,02 


38 


51 


25, 


,4 


-24.403 


2 


.273 


17, 


.649 


1, 


.160 


-1 


.505 


1 


.957 


n 


91 


BOO 24 


7 


38 


18, 


,18 


38 


52 


15, 


,9 


-21.680 


1 


.110 


17, 


.544 


1, 


.158 


-3 


.130 


1 


.991 


n 


92 


mask 2 


7 


38 


15, 


,58 


38 


51 


28, 


,4 


-15.183 


2 


.305 


18, 


.725 


0, 


.977 


-1 


.917 


1 


.993 


n 


93 


mask 1 


7 


38 


15, 


,59 


38 


54 


22, 


,0 


-18.895 


2 


.293 


18, 


.275 


1, 


.083 


-1 


.837 


2 


.001 


n 


94 


mask 2 


7 


38 


17 


,56 


38 


51 


57 


,5 


-18.071 


2 


.288 


18, 


.286 


1, 


.036 


-1 


.815 


2 


.004 


n 


95 


mask 1 


7 


38 


18, 


.48 


38 


53 


32, 


.0 


-20.034 


2 


.577 


19, 


.910 


0, 


.968 


-1 


.820 


2, 


.036 


y 


95 


mask 2 


7 


38 


18, 


,48 


38 


53 


32, 


.0 


-15.769 


2 


.452 


19, 


.910 


0, 


.968 


-1, 


.660 


2 


.036 


y 


96 


mask 1 


7 


38 


18, 


.90 


38 


52 


31, 


.8 


-18.255 


2 


.267 


17 


.237 


1, 


.204 


-1, 


.885 


2 


.058 


n 


97 


mask 2 


7 


38 


12, 


.63 


38 


51 


01, 


.2 


48.942 


2 


.719 


19, 


.990 


0, 


.814 


-0 


.742 


2 


.058 


n 


98 


mask 2 


7 


38 


17, 


.64 


38 


51 


50, 


.6 


-17.449 


2 


.334 


19 


.055 


1, 


.020 


-1, 


.735 


2 


.075 


n 


99 


mask 1 


7 


37 


58, 


,86 


38 


53 


56, 


,7 


-17.789 


2 


.334 


17, 


.790 


1, 


.077 


-1 


.182 


2 


.142 


n 


100 


mask 1 


7 


38 


17 


,81 


38 


54 


06, 


,1 


46.917 


2 


.283 


17, 


.032 


1, 


.201 


-1 


.329 


2 


.164 


n 


101 


mask 2 


7 


37 


57 


.76 


38 


53 


29, 


,3 


-16.325 


2 


.269 


16, 


.959 


1, 


.202 


-1 


.473 


2 


.169 


n 


102 


BOO III89 


7 


38 


09, 


.79 


38 


50 


45, 


,4 


-26.050 





.920 


15, 


.761 


1, 


.507 


-3 


.515 


2 


.173 


n 


103 


B09 III107 


7 


38 


15, 


,55 


38 


51 


13, 


,5 


-17.430 





.990 


17, 


.156 


1, 


.230 


-3 


.215 


2 


.176 


n 


104 


BOO 12 


7 


38 


18, 


,63 


38 


53 


52, 


,0 


-17.700 


1 


.160 


17, 


.409 


1, 


.186 


-3 


.159 


2 


.187 


y 


104 


mask 2 


7 


38 


18, 


,63 


38 


53 


52, 


,1 


-19.898 


2 


.282 


17, 


.409 


1, 


.186 


-1 


.679 


2 


.188 


y 


105 


mask 1 


7 


37 


57 


,27 


38 


52 


58, 


,0 


-24.719 


2 


.331 


17 


.816 


1, 


.122 


-1 


.915 


2 


.188 


y 


105 


mask 2 


7 


37 


57 


,27 


38 


52 


58, 


.0 


-25.534 


2 


.372 


17, 


.816 


1, 


.122 


-1 


.864 


2 


.188 


y 


106 


mask 1 


7 


38 


19, 


.64 


38 


53 


15, 


.6 


-28.739 


2 


.271 


15, 


.909 


1, 


.471 


-1 


.774 


2 


.193 


n 


107 


mask 1 


7 


38 


18, 


.16 


38 


54 


04, 


.2 


-16.560 


2, 


.557 


19, 


.591 


0, 


.982 


-1, 


.633 


2, 


.205 


n 


108 


mask 1 


7 


38 


16, 


.77 


38 


51 


22, 


.6 


-17.489 


2 


.421 


19, 


.347 


1, 


.006 


-1, 


.502 


2 


.225 


y 


108 


mask 2 


7 


38 


16, 


.77 


38 


51 


22, 


.6 


-20.194 


2 


.407 


19, 


.347 


1, 


.006 


-1, 


.481 


2 


.225 


y 


109 


B09 11118 


7 


38 


16, 


.93 


38 


54 


25, 


.2 


-20.750 


1, 


.040 


16, 


.849 


1, 


.242 


-3 


.294 


2 


.225 


n 


110 


mask 1 


7 


38 


20, 


.01 


38 


52 


48, 


.8 


-26.997 


2 


.290 


17 


.987 


1, 


.103 


-1, 


.522 


2 


.240 


n 


111 


mask 1 


7 


38 


17 


.13 


38 


51 


24, 


.5 


-19.471 


2 


.383 


19, 


.084 


1, 


.021 


-1 


.704 


2 


.255 


n 


112 


mask 1 


7 


37 


57 


,88 


38 


51 


56, 


,9 


-18.050 


2 


.280 


16, 


.595 


1, 


.305 


-1 


.596 


2 


.284 


y 


112 


mask 2 


7 


37 


57 


,88 


38 


51 


56, 


,9 


-17.875 


2 


.480 


16, 


.595 


1, 


.305 


-1 


.570 


2 


.284 


y 


113 


mask 1 


7 


37 


57 


,57 


38 


51 


58, 


,9 


-20.987 


2 


.392 


19, 


.232 


1, 


.017 


-1 


.865 


2 


.325 


n 


114 


mask 1 


7 


38 


20, 


,54 


38 


53 


06, 


,4 


-12.780 


2 


.510 


19, 


.812 


0, 


.959 


-1 


.755 


2 


.349 


n 


115 


mask 2 


7 


38 


20, 


,23 


38 


52 


12, 


,1 


-19.165 


2 


.285 


17, 


.866 


1, 


.130 


-1 


.776 


2 


.390 


n 


116 


B09 III86 


7 


38 


09, 


,54 


38 


50 


31, 


,7 


-18.890 





.990 


16, 


.973 


1, 


.248 


-3 


.259 


2 


.395 


n 


117 


mask 1 


7 


37 


58, 


,19 


38 


51 


36, 


,5 


-16.915 


2 


.710 


19 


.935 


0, 


.961 


-1 


.725 


2 


.396 


n 


118 


mask 2 


7 


38 


18, 


.02 


38 


51 


23, 


.1 


-19.906 


2 


.279 


18, 


.138 


1, 


.103 


-1 


.901 


2 


.402 


n 


119 


mask 1 


7 


37 


56, 


.40 


38 


53 


28, 


.0 


-15.372 


2 


.296 


17 


.086 


1, 


.238 


-1, 


.725 


2 


.420 


y 


119 


mask 2 


7 


37 


56, 


.40 


38 


53 


28, 


.0 


-15.207 


2 


.271 


17 


.086 


1, 


.238 


-1, 


.726 


2 


.420 


y 


120 


mask 1 


7 


38 


15, 


.28 


38 


54 


57, 


.7 


-24.280 


2 


.747 


19, 


.578 


0, 


.982 


-1, 


.197 


2 


.434 


n 


121 


mask 2 


7 


37 


56, 


.68 


38 


53 


54, 


.3 


-20.776 


2 


.284 


17 


.126 


1, 


.202 


-1 


.794 


2 


.506 


n 


122 


mask 1 


7 


38 


20, 


,92 


38 


53 


44, 


,1 


-21.836 


2 


.331 


18, 


.922 


1, 


.042 


-1 


.851 


2 


.550 


n 


123 


mask 1 


7 


37 


55, 


,40 


38 


52 


46, 


,0 


-18.686 


2 


.436 


19, 


.246 


1, 


.037 


-1 


.682 


2 


.556 


n 


124 


mask 2 


7 


37 


56, 


,48 


38 


51 


53, 


,3 


-18.725 


2 


.276 


17, 


.698 


1, 


.160 


-1 


.797 


2 


.557 


n 


125 


mask 1 


7 


37 


55, 


,80 


38 


52 


13, 


,5 


-21.527 


2 


.295 


17, 


.831 


1, 


.145 


-1 


.753 


2 


.568 


y 


125 


mask 2 


7 


37 


55, 


,80 


38 


52 


13, 


,5 


-19.102 


2 


.291 


17, 


.831 


1, 


.145 


-1 


.644 


2 


.568 


y 


126 


mask 2 


7 


38 


01, 


,20 


38 


50 


38, 


,9 


-21.521 


2 


.287 


17, 


.898 


1, 


.138 


-1 


.906 


2 


.676 


n 


127 


mask 1 


7 


37 


54, 


,67 


38 


53 


00, 


,3 


-20.692 


2 


.506 


19, 


.680 


1, 


.000 


-1 


.711 


2 


.695 


n 


128 


mask 2 


7 


38 


08, 


,56 


38 


50 


13, 


,0 


-18.569 


2 


.453 


19 


.612 


0, 


,901 


-1 


.388 


2 


.698 


n 


129 


B09 II59 


7 


38 


13, 


.97 


38 


55 


26, 


.7 


-23.200 


1 


.110 


17, 


.641 


1, 


.164 


-3 


.102 


2, 


.744 


n 


130 


B09 IV93 


7 


37 


56, 


.45 


38 


51 


22, 


.8 


-20.640 





.690 


16, 


.456 


1, 


.341 


-3 


.373 


2 


.805 


n 


131 


mask 2 


7 


38 


22, 


.80 


38 


52 


28, 


.1 


-30.504 


2 


.269 


17 


.466 


1, 


.177 


-1, 


.786 


2 


.817 


n 


132 


mask 2 


7 


38 


21, 


.96 


38 


51 


46, 


.7 


-23.469 


2 


.414 


19, 


.638 


0, 


.971 


-1, 


.733 


2 


.854 


n 


133 


mask 1 


7 


38 


23, 


.62 


38 


52 


30, 


.7 


-21.792 


2 


.675 


19 


.822 


0, 


.967 


-1 


.652 


2 


.968 


y 


133 


mask 2 


7 


38 


23, 


,62 


38 


52 


30, 


,7 


-22.507 


2 


.615 


19 


.822 


0, 


,967 


-1 


.320 


2 


.968 


y 


134 


mask 1 


7 


38 


19, 


,20 


38 


50 


47, 


,8 


-19.527 


2 


.270 


17 


.365 


1, 


,204 


-1 


.839 


2 


.969 


y 


134 


mask 2 


7 


38 


19, 


,20 


38 


50 


47, 


.8 


-19.462 


2 


.266 


17, 


.365 


1, 


.204 


-1 


.814 


2 


.969 


y 


135 


mask 1 


7 


37 


54, 


,00 


38 


54 


02, 


.5 


-18.755 


2 


.347 


18, 


.664 


1, 


.057 


-1 


.833 


3 


.040 


n 


136 


B09 III37 


7 


38 


13, 


,72 


38 


50 


00, 


,5 


-21.300 


1 


.200 


17, 


.092 


1, 


.200 


-3 


.240 


3 


.079 


n 


137 


mask 2 


7 


38 


14, 


,52 


38 


50 


00, 


,2 


-4.134 


2 


.292 


17, 


.985 


1, 


.071 


-1 


.160 


3 


.138 


n 


138 


mask 2 


7 


38 


22, 


,74 


38 


51 


26, 


,0 


-20.773 


2 


.298 


18, 


.517 


0, 


.987 


-1 


.935 


3 


.141 


n 


139 


B09 148 


7 


37 


58, 


,06 


38 


55 


22, 


,2 


-18.470 





.710 


16, 


.911 


1, 


.220 


-3 


.283 


3 


.188 


n 


140 


B09 1105 


7 


37 


54, 


.13 


38 


54 


32, 


.7 


-21.430 





.970 


16, 


.892 


1, 


.255 


-3 


.279 


3 


.238 


n 
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141 


mask 2 


7 


38 


08 


.76 


38 


49 


40 


.2 


-21.026 


2, 


.503 


19. 


.655 





.980 


-1, 


.312 


3.245 


n 


142 


mask 2 


7 


38 


18 


.47 


38 


50 


18 


.3 


-19.292 


2, 


.560 


19. 


.985 





.930 


-1, 


.349 


3.251 


n 


143 


mask 1 


7 


37 


51, 


.26 


38 


53 


07 


.4 


-22.107 


2, 


.317 


18. 


.272 


1, 


.122 


-1, 


.837 


3.363 


n 


144 


mask 1 


7 


38 


23 


.06 


38 


51 


04 


.5 


-20.391 


2, 


.322 


18. 


.527 


1, 


.033 


-1, 


.551 


3.377 


y 


144 


mask 2 


7 


38 


23 


.06 


38 


51 


04 


.5 


-22.086 


2, 


.316 


18, 


,527 


1 


.033 


-1, 


.463 


3.377 


y 


145 


mask 1 


7 


38 


26 


.02 


38 


52 


42 


.1 


-20.035 


2, 


.351 


18, 


,872 


1 


.016 


-1, 


.596 


3.414 


n 


146 


mask 2 


7 


38 


18 


.19 


38 


49 


52 


.9 


-21.514 


2, 


.314 


18, 


,866 


1 


.021 


-1, 


.591 


3.571 


n 


147 


B09 II23 


7 


38 


09 


.58 


38 


56 


30 


.0 


-20.670 


0, 


.890 


16, 


,655 


1 


.286 


-3, 


.334 


3.591 


n 


148 


mask 


2 


7 


38 


23 


.89 


38 


50 


51 


.3 


-22.526 


2, 


.270 


17, 


,190 


1 


.211 


-1, 


.871 


3.634 


n 


149 


mask 


1 


7 


38 


24 


.59 


38 


54 


50 


.4 


-21.273 


2, 


.468 


19, 


,549 





.989 


-1, 


.666 


3.673 


n 


150 


mask 


1 


7 


37 


53 


.50 


38 


50 


35 


.5 


-31.810 


2, 


.828 


19, 


,942 





.959 


-1, 


.598 


3.733 


n 


151 


mask 


2 


7 


38 


26 


.89 


38 


51 


30 


.7 


-23.012 


2, 


.380 


19, 


,129 


1 


.020 


-1, 


.713 


3.843 


n 


152 


mask 


1 


7 


38 


25, 


.26 


38 


55 


15, 


.9 


16.419 


2, 


.626 


19, 


,655 


0, 


.978 


-0, 


.873 


4.018 


n 


153 


mask 


1 


7 


37 


48, 


.61 


38 


51 


20 


.6 


-22.418 


2, 


.325 


18. 


.673 





.999 


-1, 


.610 


4.180 


n 


154 


mask 


2 


7 


38 


27 


.15 


38 


50 


43, 


.3 


-18.856 


2, 


.522 


19. 


.915 





.962 


-1, 


.540 


4.240 


n 


155 


mask 


1 


7 


37 


46 


.06 


38 


51 


41, 


.4 


-21.155 


2, 


.281 


17. 


.821 


1, 


.094 


-1, 


.617 


4.538 


n 


156 


mask 


2 


7 


38 


19 


.92 


38 


48 


50 


.4 


18.899 


2, 


.518 


19, 


,661 





.867 


-0, 


.709 


4.641 


n 


157 


mask 


1 


7 


37 


45 


.44 


38 


54 


23 


.6 


-20.425 


2, 


.356 


18, 


,950 


1 


.015 


-1, 


.659 


4.726 


n 


158 


mask 


1 


7 


38 


31 


.48 


38 


51 


12 


.1 


-21.597 


2, 


.277 


17, 


,598 


1 


.172 


-1, 


.772 


4.788 


n 


159 


mask 


1 


7 


38 


33 


.01 


38 


53 


45 


.2 


-21.967 


2, 


.264 


16, 


,379 


1 


.304 


-1, 


.452 


4.841 


n 


160 


mask 


1 


7 


38 


32 


.27 


38 


54 


43 


.4 


-20.964 


2, 


.292 


17, 


,018 


1 


.233 


-1, 


.530 


4.964 


n 


161 


mask 


1 


7 


37 


42 


.49 


38 


51 


11 


.9 


-14.861 


2, 


.338 


19, 


,085 


1 


.020 


-1, 


.769 


5.348 


n 


162 


mask 


1 


7 


37 


40 


.59 


38 


52 


03 


.0 


13.429 


2, 


.273 


16, 


,840 


1 


.190 


-1, 


.277 


5.502 


n 


163 


mask 


2 


7 


38 


31 


.98 


38 


49 


47 


.7 


-21.679 


2, 


.429 


19, 


,249 


1 


.006 


-1, 


.545 


5.533 


n 


164 


mask 


1 


7 


37 


40 


.96 


38 


51 


05 


.7 


-21.438 


2, 


.260 


16, 


,637 


1 


.285 


-1, 


.909 


5.663 


n 


165 


mask 


2 


7 


38 


32, 


.03 


38 


49 


13, 


.9 


-52.154 


2, 


.368 


18, 


,451 


1, 


.100 


-0, 


.993 


5.877 


n 


166 


mask 


1 


7 


38 


41, 


.05 


38 


52 


26 


.0 


-18.550 


2, 


.322 


17. 


.867 


1, 


.053 


-2, 


.047 


6.351 


n 


167 


mask 


2 


7 


38 


30 


.47 


38 


47 


46, 


.5 


-21.544 


2, 


.512 


19. 


.722 





.947 


-1, 


.514 


6.686 


n 


168 


mask 


1 


7 


37 


33 


.01 


38 


51 


39 


.9 


-25.755 


2, 


.603 


19. 


.637 





.950 


-1, 


.705 


7.022 


n 


169 


mask 


2 


7 


38 


42, 


.77 


38 


49 


41, 


.0 


-160.538 


2, 


.735 


19. 


.925 





.828 


-0, 


.867 


7.412 


n 


170 


mask 


1 


7 


37 


29 


.30 


38 


52 


28 


.7 


-19.041 


2, 


.286 


17, 


,645 


1 


.119 


-1, 


.941 


7.644 


n 


171 


mask 


2 


7 


38 


34 


.78 


38 


47 


04 


.1 


-21.921 


2, 


.525 


19, 


,273 





.956 


-0, 


.822 


7.769 


n 


172 


mask 


2 


7 


38 


48 


.16 


38 


48 


19 


.6 


-5.179 


2, 


.499 


19, 


,690 





.968 


-1, 


.272 


8.981 


n 


173 


mask 


1 


7 


37 


22 


.10 


38 


52 


45 


.5 


-32.019 


2, 


.356 


18, 


,114 


1 


.155 


-0, 


.977 


9.033 


n 


174 


mask 


1 


7 


37 


22 


.89 


38 


51 


05 


.9 


-14.079 


2, 


.506 


18, 


,711 


1 


.120 


-1, 


.519 


9.064 


n 


175 


mask 


2 


7 


38 


43 


.77 


38 


46 


01 


.5 


-50.217 


2, 


.278 


17, 


,046 


1 


.239 


-1, 


.127 


9.728 


n 


176 


mask 


2 


7 


38 


49 


.47 


38 


46 


22 


.9 


-49.256 


2, 


.285 


17, 


,281 


1 


.238 


-1, 


.082 


10.311 


n 


177 


mask 


2 


7 


38 


55 


.49 


38 


47 


44 


.3 


-21.062 


2, 


.540 


16, 


,729 


1 


.317 


-1, 


.052 


10.511 


n 


178 


mask 


2 


7 


38 


56, 


.74 


38 


43 


54, 


.3 


-22.239 


2, 


.533 


19, 


,340 


0, 


.876 


-1, 


.066 


13.018 


n 



